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SECTION  I 


SUMMARY 


A.  INTRODUCTION 


Project  VT/4053  ie  directed  toward  continued  deve’opment 
of  array  processing  technology  for  nuclear  surveillance  and  exploitation  of 
the  superior  data  available  from  arrays  for  analysis  of  distant  P-waves. 

Work  during  the  period  covered  by  this  report  has  included: 

•  Processing  of  vertical  array  data  for  the  purposes  of 
evaluating  MCF  systems  and  decomposing  the  ambient 
noise  field  into  constituent  modes 

•  Theoretical  investigation  of  the  use  of  horizontal 
seismometers  in  combination  with  vertical  seismometers 
for  extraction  of  signals 

•  Processing  of  CPO  teleseisms  for  the  purposes  of  investi¬ 
gating  crustal  reverberation  deconvolution,  determining 
depth  of  focu!  stud',  li._  CPO  crustal  structure,  and 
investigating  propagation  mechanisms 

e  Further  study  of  CPO  partial  arra>  with  particular 
emphasis  on  the  problems  associated  with  seismometer 
inequali  nation 

•  A  study  of  non-directional  array  processing  at  TFO 
utilizing  the  large  cross  array, 

•  A  study  of  seismic  signal  detection  as  a  separate  problem 
from  the  problem  of  signal  extraction 

e  A  study  of  the  ability  of  the  Montana  LASA  to  seoarate 
nuclear  blast  signals  from  natural  earthquakes,  and  a 
study  of  an  operational  LASA  signal  processing  system 

•  Collection  and  preparation  of  seismic  data  to  be  used  in 
future  experimental  work 

•  Development  of  an  automated  mapping  program  which 
makes  it  practical  to  present  the  world  in  wavenumber 
space  relative  to  a  particular  station 

•  Investigation  of  the  predictability  of  seismic  noise  in  the 
shallow-buried  UBO  array 

•  A  study  into  the  extraction  of  teleseiamic  signals  from 
cultural  noise  af  UBO 
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B.  PROCESSING  OF  VERTICAL  ARRAY  DATA 


Section  II  of  this  report  presents  results  obtained  from  pro¬ 
cessing  vertical  array  data.  Those  results  considered  most  pertinent  to  the 
nuclear -detection  problem  are  presented  very  briefly  in  the  following  para¬ 
graphs. 

Analysis  of  vertical  array  data  has  lent  support  to  the  hypothe¬ 
sis  tnat  h’gh-velocity  body-wave  noise  is  a  significant  fraction  of  the  ambient 
noise  at  quiet  sites  in  the  frequency  range  from  0.  5  cps  1.2  cps.  An  exam¬ 
ple  is  the  modal  analysis  of  ambient  noise  at  UBO,  illustrated  in  Figure  I-i. 
The  vertical  array  analysis  at  UBO  has  also  answered  the  question  as  to  whe¬ 
ther  the  high-velocity  noise  might  originate  at  the  surface  in  the  vicinity  of 
the  array.  At  UBO,  that  20  percent  of  the  ambient  noise  power  at  1.  0  cps 
which  was  observed  on  wavenumber  analysis  of  the  surface  array  to  be  high- 
velocity  noise,  has  now  been  identified  on  i.he  vertical  array  as  body  waves 
originating  below  the  vertical  array.  Further  analysis  at  UBO  and  a  granite 
well  have  supported  the  hypothesis  that  the  body-wave  noise  component  is 
relatively  uniform  from  one  site  to  the  next. 

This  hypothesis  explains  the  poor  signal-to-noise  improvement 
obtained  in  this  frequency  range  at  very  quiet  sites  by  limited  aperture  sur¬ 
face  arraye  and  by  deeply  buried  seismometers  (down  boreholes).  It  al;  o 
limits  the  steady  state  performance  of  vertical  arrays  at  the  UBO  and  granite 
well  sites  to  a  fev,'  db  in  this  frequency  range.  Such  systems  must  be  consid¬ 
ered  to  be  limited  to  reduction  of  trapped-mode  noise  bursts  at  very  quiet 
sites  in  this  frequency  range.  Performance  for  such  noise  bursts  should  be 
related  to  the  theoretical  analyses  which  have  been  previously  published. 
Evaluation  of  vertical  array  performance  in  this  sense  has  not  been  completed. 

In  the  frequency  range  below  0,  3  cps,  most  evidence  supports 
the  hypothesis  that  Rayleigh- wave  noise  is  dominant.  Although  the  vertical 
arrays  examined  have  limited  aperture  for  this  very  low  frequency  noise, 
theoretical  studies  suggested  a  possible  moderate  performance  at  very  low 
frequencies.  Steady  state  noise  reduction  cf  5  to  6  db  has  been  observed  at 
the  very  low  frequencies  at  the  granite  well  site.  However,  this  could  be  due 
to  instrument  gain  problems  and  requires  further  examination. 

Analysis  suggests  (Figure  1-1}  that  above  1.0  to  1.2  cps  the 
dominant  noise,  even  at  quiet  sites,  is  dominantly  trapped-mode  noise  or  noise 
which  is  not  sufficiently  correlated  at  the  experimental  ipacings  to  be  identi¬ 
fiable.  In  the  1.  5  cps  to  3.  C  epo  band,  gains  of  JW  or  better  have  been 
obtained  at  quiet  sites  such  as  TFO  and  UBO  and  a.  noisy  sites  such  as  CPO, 
with  small  aperture  surface  arrays  of  N  elements.  In  this  frequency  range, 
it  is  reasonable  to  expect  good  signal-to-noise  improvement  from  vertical 
arrays  with  the  limitations  suggested  by  previous  theoretical  work.  Improve¬ 
ments  of  the  order  of  10  db  have  been  found  at  UBO  in  this  frequency  range. 
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Significant  lesser  gains  have  been  observed  at  the  granite  well,  although 
data  analyzed  to  date  has  probably  been,  contaminated  by  tape  noise  in  the 
high-frequency  rang<?t 

In  evaluation  of  vertical  arrays  in  the  frequency  range  below 
0.  5  cps  and  above  1*  5  cps,  results  to  date  have  been  less  than  conclusive 
because  of  difficulties  in  the  specification  of  the  signal  model.  Use  of 
"measured  signal"  models  is  limited  to  the  frequency  band  where  good  tele- 
seism  signal-to-noise  ratios  have  been  obtained  on  individual  seismometers, 
which  thus  fax  have  not  extended  to  the  low  and  high  frequencies  of  interest. 
Theoretical  models  are  harassed  by  uncertainties  in  the  relative  gains  of  the 
different  seismometers.  The  best  approach  would  appear  to  be  the  use  of  a 
theoretical  signal  model  with  limited  gain  slop  ,  using  the  beat  available 
estimate  of  relative  signal  gain.  This  approach  has  not  yet  been  adequately 
tested  at  any  of  the  vertical  array  sites. 

The  major  potential  value  of  vertical  arrays  in  VELA  UNIFORM 

is  twofold: 

(1)  Reduction  of  noise  in  the  0.  5  to  1.  5  cps  range  to  body  wave  noise. 

At  quiet  sites,  this  imporvement  will  be  small  on  an  average  basis. 
However,  noise  "bursts'^  or  intermittent  "local"  noise  should  be 
reduced. 

(?,}  The  most  significant  potential  cain  lies  in  the  frequency  band  above 
1.  5  cps.  It  may  turn  out  that  smaller  spacings  than  those  employed 
so  far  may  be  required  for  maximum  vertical  array  performance  in 
this  high  frequency  range. 

The  principal  direction  of  the  continuing  vertical  array  work 
is  in  the  evaluation  of  performance  above  1.5  cps.  If  performance  signifi¬ 
cantly  better  than  the  observed  10  db  at  UBO  can  be  achieved  at  quiet  sites, 
a  possibly  attractive  system  may  be  a  moderate  depth  vertical  array  for  high- 
frequency  performance  that  is  combined  with  some  horizontally-displaced, 
near-surface  radial  seismometers  to  handle  the  low-frequency  microseism 
peak.  The  latter  possibility  is  discussed  in  the  following  section. 


C,  SYSTEMS  OF  HORIZONTAL  AND  VERTICAL  SEISMOMETERS 

Special  Report  No.  7  entitled  "Array  Research  Theoreti¬ 
cal  Capability  of  Systems  of  Horizontal  Seismometers  for  Predicting  a 
Vertical  Component  in  Ambient  Trapped  Mode  Noise"  was  published 
9  November  1965.  A  discussion  of  that  report  follows. 
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i.  Multicomponent  Seismometers 


Mathematical  formulae  were  derived  for  the  noise-rejection 
properties  of  a  2-component  seismometer,  consisting  of  1  vertical-com¬ 
ponent  and  1  horizontal-component  instrument  at  the  same  point.  For 
most  applications,  this  system  may  be  considered  equivalent  to  a 
3-component  seismometer.  It  is  possible  to  design  a  processing  system 
wherein  different  filters  are  applied  to  the  outputs  of  2  horizontals  and  the 
noise  estimate  is  obtained  by  summing  the  filter  outputs. 


Some  general  conclusions  were  stated  regarding  the 
usefulness  of  a  single  horizontal-component  instrument  for  eliminating 
noise  from  the  output  of  a  vertical  instrument  at  the  same  location.  It 
appeared  that  such  effects  as  system  noise,  gain  fluctuations  and  Love  waves 
would  not  generally  present  serious  problems,  although  in  some  cases  they 
may  very  well  do  so.  The  most  important  consideration  is  usually  the 
properties  of  the  noise  to  be  eliminated.  If  significant  noise  rejection  is  to 
be  achieved,  then  some  rather  stringent  conditions  must  be  satisfied.  Even 
when  therre  is  only  a  single-noise  mode  present,  it  is  essential  that  the 
azimuths  of  the  noise  sources  be  confined  within  a  range  somewhat  less 
than  180°  and  that  the  horizontal  be  oriented  near  the  center  of  this  range. 
This  requirement  follows  from  the  obvious  fact  that  a  filter  which  cancels 
out  noise  traveling  in  a  given  direction  must  also  amplify  noise  traveling 
in  the  opposite  direction.  Ip.  practice,  this  requirement  may  be  partially 
met  in  major  microseismic  "storms." 

If  more  than  one  noise  mode  is  important,  large  differences 
in  the  horizontal-to-vertical  transfer  functions,  k(t),  cannot  be  permitted. 

A  two-mode  situation  was  presented  in  which  no  noise  rejection  is  possible 
because  one  mode  is  characterized  by  a  progade  particle  motion  and  the 
other  is  associated  with  retrograde  motion. 

With  a  3 -component  seismometer  at  point  location,  it  is 
theoretically  possible  to  determine  both  the  direction  of  propagation  and  the 
shape  of  the  particle  orbit  for  an  observed  Rayleigh  wave.  However,  the 
sense  of  propagation  and  the  Bense  of  the  orbital  motion  cannot  be  found 
unless  a  spatial  separation  is  introduced  into  the  system.  For  this  reason 
it  is  concluded  that  ordinary  3 -component  systems  offer  little  promise  for 
most  applications  and  that  attention  must  be  directed  toward  arrays  A 
horizontal-  and  vertical-component  seismometers. 
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2.  Multicomponent  Arrays 


Theoretical  results  obtained  in  this  study  indicated  that  arrays 
of  horizontal -component  seismometers  should  prove  to  be  useful  tools  for 
the  removal  of  trapped-mode  noise  from  the  outputs  of  vertical -component 
instruments.  A  necessary  condition  is  that  the  horizontals  must  be 
spatially  separated  from  the  vertical.  In  the  case  of  a  single-noise  mode, 
arrays  with  a  single  ring  of  horizontal  seismometers  should  perform  best 
when  the  array  diameter  is  approximately  one -half  wavelength.  Signal 
enhancement  systems  employing  only  vertical  instruments  require,  in 
general,  array  diameters  of  at  least  one  wavelength.  Thus,  multicompon¬ 
ent  arrays  should  offer  meaningful  advantages  in  terms  of  land  and 
telemetry  requirements. 


It  is  unlikely  that  difficulties  presented  by  system  noise, 
uncorrelated  seismic  noise  and  Love  waves  should  be  any  more  serious 
than  they  are  in  the  case  of  vertical-component  arrays.  The  effect  of 
additional  noise  modes  on  multicomponent  array  performance  has  not  been 
studied  yet.  However,  it  is  reasonable  to  assume  that  a  multiplicity  of 
modes  may  be  dealt  with  by  the  application  of  multichannel  filter  techniques 
to  the  outputs  of  multiple  rings  of  Horizontals.  It  has  been  shown  that  an 
array  consisting  of  a  number  of  rings  of  verticals  can  be  useful  in  the 
presence  of  a  similar  number  of  noise  modes,  and  there  is  no  known  reason 
for  assuming  that  this  usefulness  might  be  a  property  peculiar  to  verlical- 
component  arrays. 

The  results  derived  theoretically  for  multicomponent 
seismometer  arrays  are  quite  conducive  to  optimism.  It  was  shown  that, 
if  the  array  dimensions  are  suitable  (r/X~ 0.25),  the  theoretical  prediction 
error  is  less  than  0.01,  the  Love  wave  response  function  is  less  than 
0,01  and  the  Rayleigh  wave  response  function  is  approximately  0.7.  Thus, 
it  is  theoretically  possible  to  predict  (and  hence  remove)  more  than 
99  percent  of  the  Rayleigh  noise  on  the  vertical  component  in  the  frequency 
range  in  which  the  wavelength  is  appropriate.  The  low  value  of  the  Love 
wave  response  function  implies  that  very  littlo  extraneous  noise  power 
should  be  introduced  into  the  system  output  as  a.  result  of  Love  waves 
appearing  in  the  outputs  of  the  horizontal  instruments.  The  high  value  of 
the  Rayleigh  wave  response  function  implies  that  the  filter  responses 
need  not  be  unduly  large  and  that  uncorrelated  noise  from  the  horizontals 
will  not  be  amplified  to  a  serious  degree.  The  finding  that  ail  three  response 
functions  take  on  desirable  values  in  the  same  range  of  wavelength  is  both 
unexpected  and  fortuitous. 


1-6 


From  the  results  of  the  theoretical  investigations  presented 
in  this  report,  it  is  concluded  that  multicomponent  seismometer  arrays 
offer  a  great  deal  of  promise  for  signal  enhancement  applications.  Experi¬ 
mental  investigations  of  such  systems  is  underway  using  data  from  1  FO 
It  is  recommended  that  experimental  investigations  of  such  systems  be 
undertaker,  at  a  site  where  a  largo  percentage  of  the  ambient  noise  is 
trapped-mode  noise. 

D.  ANALYSIS  OF  CPO  TELESEISMS 

Several  large  ensembles  of  CPO  teleseis.ns  have  been 
analyzed  for  the  purposes  of: 

•  Investigating  crustal  reverberation  deconvolution 

•  Determining  depth  of  focus 

•  Studying  CPO  crustal  structure 

•  Investigating  propagation  mechamisms 

0  Investigating  teleseismic  signal  properties 
relevant  to  the  nuclear  detection  problem 

Some  of  these  investigations  have  been  discussed  in  Semi¬ 
annual  Technical  Reports  2  and  3.  Some  further  discussion,  particularly 
with  regard  to  CPO  crustal  structure,  if  provided  in  Section  HI.  A  com¬ 
plete  report  on  these  studies  will  be  provide  in  Special  Report  14.  Fairly 
good  evidence  of  a  significant  shallow  layering  effect  on  teleseisms  at  CPO 
has  been  found.  Compensation  for  tills  effect  by  deconvolution  has  given 
improvement  in  some  senses  and  is  probably  not  highly  range -sensitive.  The 
deeper  reflectors,  the  Conrad  and  Mohorovicic,  do  not  present  a  significant 
signal  alteration  oroblem  in  the  short  period  range  and  have  been  only 
marginally  identified  thus  far  in  teleseism  autocorrelations.  There  is  some 
question  as  to  the  adequacy  of  the  existing  ensemble  statistics  to  unambigu¬ 
ously  identify  deep  reflectors  due,  in  part,  to  the  range  and  azimuth  depend¬ 
ence  of  the  delay  time  for  the  deep  reflectors.  It  is  of  interest  that  the 
crustal  effects  are  clearly  most  evident  on  deep  focus  earthquakes  a?  a  res  lit 
of  the  source  simplicity  of  deep-focus  events  as  compared  with  shallow  events. 
It  is  not  unlikely  that,  even  in  view  of  the  capabilities  of  ensemble  statistics, 
studies  of  station  response  and  propagation  path  phenomena  should  be  more 
fruitful  when  based  on  a  small  ensemble  of  deep-focus  events  rather  tnan  a 
large  ensemble  of  shallow  events. 

E.  FURTHER  STUDY  OF  CPO  PARTIAL  ARRAYS 

The  seismometer  equalization  problem  at  CPO  was  considered 
and  the  results  indicated  that  conclusions  previously  drawn,  regarding  the 
effectiveness  of  partial  arrays  at  CPO,  remain  valid. 
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F.  nondirec tional  cross  array  processing  at  tfo 

An  MCF  system  designed  to  operate  on  the  large  cross  array 
at  TFO  provided  3  db  to  5  db  greater  S/N  improvement  than  a  simple  summa¬ 
tion  at  frequencies  below  about  1.  0  cps.  At  0.2  cps,  the  MCT  system  Pro" 
vided  8-db  S/N  improvement  relative  to  a  simple  summation.  Above  about 
1.  0  cp  ,  the  MCF  and  summation  provided  about  equal  S/N  improvement. 

Results  obtained  are  consistent  with  the  hypothesis  that 
trapped-mode  noise  at  TFO  is  3  to  4  db  above  body-wave  noise  in  the  0.25 
to  0.  75  cps  range  ana  about  13  db  above  body-wave  noise  in  the  0,  15  to  0.25 
cr  s  range.  It  is  interesting  that  the  latter  result  is  in  semiquantitative  agree¬ 
ment  with  the  coherence  analyses  reported  in  Section  VI  of  Semiannual  Tech¬ 
nical  Report  3.  The  wavenumber  response  analysis  is  of  interest  in  the 
general  LASA  subsystem  problem.  In  the  range  of  0.  15  to  0,25  cps,  the 
array  is  clearly  concentrating  on  a  trapped-mode  microseismic  wavetrain 
coming  from  N  60°  E  at  a  velocity  of  3  to  4  km/sec.  In  the  0.  5  to  1.  0  cps 
oand,  the  10-km  array  is  primarily  concentrating  on  a  high-velocity  body- 
wave  component  coming  from  the  southeast  at  a  velocity  probably  in  excess  of 
of  16  km/sec.  Above  1.  5  cps,  the  array  provides /N  improvement  or  better, 
whether  employed  as  a  simple  beam  former  or  as  an  adaptive  noise  rejector. 

G.  DETECTION  OF  SEISMIC  SIGNALS 

Special  Report  No.  8  entitled  "Array  Research  Preliminary 
Report,  Matrix  Multiply  Detection  Processing  of  Array  Data"  was  published 
13  October  1965.  Assumptions,  derivations  of  equations,  and  an  illustration 
of  the  method  of  signal  detection  by  probabilistic  processing  were  presented. 

Results  of  applying  a  probabilistic  processor  to  actual 
signals  are  presented  in  Section  VI  of  this  report.  This  processor  was 
designed  using  a  theoretical  model  of  signal  and  noise.  Three  events  were 
processed  and  the  results  were  compared  with  a  detection  scheme  employ¬ 
ing  Wiener  signal-extraction  processing  and  the  square  and  integrate 
process  on  the  output  trace.  For  one  of  the  signals,  an  indication  of 
superiority  of  *ve  matrix- multiply  technique  over  the  Wiener  technique 
was  obtained. 


It  is  planned  to  evaluate  the  matrix  multiply  detection 
technique  using  a  processor  designed  from  measured  CPO  noise  correlations 
It  is  also  planned  to  compare  the  results  with  a  comparable  Wiener  square^ 
and  integrate  process  and  with  a  technique  known  as  Analysis  of  Variance. 


*BOoker,  Aaron,  1965:  Analysis  of  variance  as  a  method  for  seismic  signal 
detection,  SDL  Rpfc.  No.  116,  Contract  AF  33(657)- 12447,  VT/Z037,  25  Feb. 


H.  LASA  PROCESSING,  DATA  COLLECTION,  AUTOMATED  MAPPING 
AND  PREDICTABILITY  OF  NOISE  FROM  SHALLOW -BURIED  ARRAYS 

A  study  was  completed  to  determine  the  effectiveneas  of 
LASA  multichannel  processing  in  separating  teleseiemic  nuclear  blast 
signals  (originating  in  particular  regions  on  the  earth)  from  teleseiemic 
earthquakes  (originating  in  certain  seismically  active  regions  of  the 
earth).  Several  different  cases  were  investigated  and  the  results  are  given 
in  Section  VII.  A  conceptual  method  of  processing  LASA  data  on-line  was 
evolved  and  is  presented  in  summary  form  in  Section  X. 

Data  collection  programs  in  support  of  Dallas -based  research 
efforts  have  been  active  and  are  summarized  in  Section  IX. 

An  automated  mapping  system  for  presenting  the  world  in 
wavenumber  space  is  presented  in  Section  XL 

The  UBO  studies  dealing  with  spatial  predictability  of 
shallow-buried  seismic  noise  and  extraction  of  signals  from  cultural  noise 
are  incomplete. 
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SECTION  II 

NEW  RESULTS  FROM  PROCESSING  OF  VERTICAL  ARRAY  DATA 

A.  RESULTS  OF  FOlSE  STUDY,  FILTER  DESIGN  AND  EVALUATION 
BASED  ON  PI2VY  DATA 

Deep -well  seismic  data  which  have  been  obtained  from  the 
granite  well  designated  PI2V/Y  are  displayed  in  Figures  II -1  and  II -2.  A3 
indicated,  noise  recordings  wnre  made  with  a  surface  instrument  and 
deep  -veil  instruments  at  four  different  depths.  Signal  recordings  were 
available  for  the  surface  instrument  and  for  two  depths.  Each  recording 
consists  of  simultaneous  data  am  the  surface  and  two  depths  m  the  well. 

Start  times  (GMT)  of  the  noise  samples  are  indicated  in 
Figure  II-l.  Noise  samples  1,  2  and  3  were  recorded  between  7  a.m.  and 
11  a.m.  (local  time)  while  noise  sample  4  was  recorded  near  midnight 
(local  time).  Differences  in  noise  power  may  indicate  a  difference  between 
day  and  night  levels,  but  also  could  result  from  a  long-term  effect  since 
the  recordings  span  an  internal  of  3—1/2  weeks.  Sources  of  the  tt ’eseismic 
signals  have  b~en  identified  as  follows: 

•  Signal  1,  Fiji  Islands  Region;  17.0°S,  177.  4°W;  origin 
time,  01/36/45.4  GMT;  C&GS  focal  depth,  386  km; 
magnitude,  5.  5  to  6.  0;  distance,  85°. 

•  Signal  2.  Aegean  Sea;  39.  4°N,  24.0°E;  origin  time, 
17/57/53.7  GMT;  C&GS  focal  depth,  18  km;  magnitude 
5.7  to  6.5;  distance,  8bs. 

Records  .  .iginally  were  digitized  with  40-  msec  sample  rate. 
After  the  application  of  a  high -cut  antialiasing  filter  (3-db  down -point  at 
3.5  cps),  the  records  were  resampled  at  120  m3ec  for  analysis.  Power 
spectra  of  the  surface  instrum  nt  outputs  (not  corrected  for  instrument 
response)  are  shown  in  Figures  II -3  ar  TI-4.  A  strong  low -frequency 
microseismic  peak  is  evident  in  all  fo.  noise  spectra,  but  does  not  appear 
at  a”  m  the  signal  spectra.  Apparently,  two  distinct  types  of  noise  have 
been  recorded.  Low-frequency  microseismic  noise  is  observed  in  short 
(less  than  10  sec)  bursts  which  seem  to  arrive  at  random  times  (Figure  II-l). 
This  type  of  noise  is  extremely  reproducible  from  one  depth  to  another  and 
its  amplitude  shows  only  slight  depth  dependence.  The  other  type  of  noise 
may  be  observed  between  the  low-frequency  bursts  and  also  before  the  signal 
onsets  (Figure  II-2)t  This  noise  has  much  more  high-frequency  energy 
and  displays  a  clear  amplitude  decrease  with  depth,  even  at  low  f.i  equencies. 
It  appears  that  this  second  type  of  noise  is  relatively  steady  and  never 
approaches  the  strength  of  the  microseismic  bursts. 
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Figure  II -2,  Teles ci3mic  Signcdi'  Recorded  in  the  Well  PI2WY 
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The  flattening  of  the  noise  spectra  above  1.  5  cps  is  an 
indication  of  system  noise  which  could  have  been  introduced  either 
in  the  recording  or  digitization  operations.  Since  system  noise  is 
dominant  in  the  high-frequency  range,  vertical  array  evaluation  is 
possible  only  for  frequencies  below  1,5  cp3. 

An  important  objective  of  this  study  was  the  evaluation  oS 
multichannel  filter  techniques  applied  to  data  from  vertical  arrays  in  wells 
such  as  PI2WY  for  the  purpose  of  signal -to -noise  ratio  improvement. 

Such  Improvement  can  be  significant  only  if  the  amplitudes  of  signal  and 
noise  depend  upon  depth  in  different  predictable  ways  or  if  the  noise  is 
not  coherent.  Deep  well -to -surface  power  spectral  ratios  computed  from 
the  available  data  are  shown  in  Figure  U-5,  Below  1,  5  cps,  only  minor 
differences  between  signal  and  noise  behavior  can  be  observed.  Coherence 
functions  computed  for  Noise  Sample  1  (Figure  11-6)  indicate  that  the  noise 
is  predictable  from  channel  to  channel  up  to  approximately  1.  0  cps  for  the 
larger  receiver  separations  (6000  ft,  8000  ft)  and  1.  5  cps  for  the  smallest 
separation  (2000  ft). 

For  multichannel  filter  design,  two  models  were  employed. 

In  both  cases,  correlation  functions  computed  from  noise  sample  1  were 
used  as  the  noise  model.  The  signal  model  for  the  first,  case  (MCF-1) 
was  obtained  by  averaging  correlation  functions  of  the  two  experimental 
signals.  For  the  second  case  (MCF-2),  a  theoretical  signal  with  white 
spectrum  in  the  range  of  0  to  3  cps  was  used.  Before  experimental 
correlations  were  computed,  each  record  was  whitened  by  applying  to  all 
three  traces  a  deconvolution  filter  designed  from  the  surface  trace.  Impulse 
responses  of  the  two  MCF  sets  are  shown  in  Figures  II -7  and  II-8. 

Results  of  applying  the  filters  to  noise  sample  1  and  signals 
1  and  2  are  shown  in  Figures  II-9  and  11-10.  For  each  filter  set,  the  output 
trace  and  the  reject  trace  (difference  between  input  trace  1  and  output  trace) 
are  shown.  For  the  high-frequency  behavior  of  the  systems  to  be  seen  more 
clearly,  a  low-cut  filter  (3-db  down-point  at  0.76  cps)  has  been  applied  to 
both  the  input  and  output  traces  (Figures  II- 11  and  11-12). 

Input  and  output  power  spectra  oi  signal  and  noise  are  shown 
in  Figure  11-13.  The  only  significant  difference  between  the  performances 
of  the  two  systems  is  that  MCF-2,  designed  on  an  assumption  of  zero- 
signal  power  above  the  frequency  of  3  cps,  rejected  everything  above  that 
frequency.  Signal-to-noise  improvements  of  the  two  systems  are  given  in 
Figure  11-14. 
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Figure  II-5a.  Deep-Well  to  Surface  Spectral  Ratio  of  Seismic  Noise 


Figure  -I -5b,  Deep -Well  to  Surface  Spectral  Ratios  of  '"eleseismic  Sisals 


TIME  (SECONDS) 
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Figure  II -1  3.  Input  and  MCF  Output  Power  Spectra  of  Noise  and  Signal 
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Figure  11-14.  S/N  Improvements  Computed  for  the  Two  MCF  Systems 


At  low  frequencies  (<  0.  5  cps)  both  systems  achieved  about 
7  to  8  db  signal -to -noise  improvement.  Since  it  is  unlikely  that  there 
could  have  been  significant  signal  energy  at  such  low  frequencies,  this 
result  implies  that  noise  in  Noise  Sample  1  was  rejected  more  effectively 
than  noise  in  the  signal  records.  Such  behavior  could  be  expected 
from  MCF-1,  designed  from  experimental  signal  statistics,  if  there 
were  a  marked  difference  in  the  character  of  the  noise.  The  failure 
of  MCF-2  to  reject  low-frequency  noise  in  the  signal  records  suggests 
that  channel -to -channel  gain  compensation  may  not  have  been  applied 
correctly,  and  that  the  apparent  changes  in  noise  character  may  be  arti¬ 
ficial  ones. 


At  intermediate  frequencies  (0.5  <  f  <  1.2  cps),  neither 
system  achieved  any  significant  improvement.  This  failure  and  the 
propeities  of  the  spectral  ratios  and  coherences  computed  for  Noise 
Sample  1  can  be  explained  most  easily  by  assuming  that  the  ambient 
noise  is  composed  mainly  of  F -waves  in  the  intermediate  frequency  range. 
Since  the  subsurface  at  the  well  site  may,  for  practical  purposes,  be 
represented  by  an  infinite  homogeneous  half  space,  the  only  kind  of  surface 
wave  which  should  exist  is  the  fundamental -mode  Rayleigh  wave.  The 
observed  peaks  and  troughs  in  the  experimental  spectral  ratios  and 
coherences  agree  very  well  with  those  found  theoretically  and  experimentally 
for  P-waves,  but  should  not  be  produced  at  all  by  Rayleigh -wave  noise. 

As  was  mentioned  above,  information  in  the  high-frequency 
range  (>  1.2  cps)  has  been  destroyed  by  system  noise.  The  lack  of 
coherence  observed  foi*  this  noise  indicates  that  it  was  probably  introduced 
in  the  analog-to-digital  conversion.  Noise  due  to  digitization  round-off 
should  be  between  30  and  40  db  lower  than  the  observed  system  noise 
level.  It  is  therefore  concluded  that  the  explanation  lies  in  inaccuracies 
in  the  analog  tape  transports.  Since  performance  in  the  intermediate- 
frequency  band  is  apparently  limited  by  the  P-wave  noise  level,  it  is 
clear  that  the  higher  frequencies  deserve  more  attention  than  they  have 
received  in  the  past.  Useful  imormation,  however,  cannot  be  obtained 
with  the  system  used  to  collect  the  data  described  here.  The  additional 
dynamic  range  afforded  by  digital  field  recording  is  sufficient  for  this 
purpose.  Alternatively,  low-cut  filtering  of  the  data  before  recording 
could  be  used  to  ease  dynamic  range  requirements. 

In  summary,  the  results  obtained  for  the  PI2WY  are 
inconclusive.  It  appears  that  useful  vertical  array  performance  must  be 
limited  to  the  frequency  range  >1.2  cps)  but  it  xs  impossible  to  estimate 
how  much,  if  any,  signal -to-noise  improvement  can  be  expected  there. 
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Table  II -1 


CATALOG  OF  SELECTED  UBO  VERTICAL  ARRAY  RECORDINGS 


Record 

Identification 

Record 

No. 

Date 

mmSSm 

Description 

Noise  Sample  1 

45 

Oct.  5 

134010 

Quiet  Noise 

Noise  Sample  2 

46 

5 

135410 

Intense  Noise 

Noise  Sample  3 

47 

5 

135810 

Quiet  Noise 

Noise  Sample  4 

50 

5 

142210 

Intermediate  -Level 
Noise 

Noise  Sample  5 

52 

5 

144410 

Intense  Noise 

Noise  Sample  6 

53 

5 

144810 

Intense  Noise 

Local  Event  1 

501 

3 

220110 

Unidentified  Local  or 
Near -Regional  Event 

Local  Event  2 

506 

5 

0Q141C 

Unidentified  Local  or 
Near -Regional  Event 

Local  Event  3 

512 

6 

002010 

Unidentified  Local  or 
Near -Regional  Event 

Signal  1 

407 

5 

002210 

Yukon  Earthquake 

Signal  2 

423 

5 

100330 

Mid -Indian  Rise 
Earthquake 

Signal  3 

403 

4 

041550 

Oregon  Offshore 
Earthquake 

Signal  * 

406 

5 

001110 

Mexico  Earthquake 

Signal  5 

412 

5 

081910 

Rat  Island  Earthquake 

Signal  6 

420 

5 

035100 

Chile  Offshore  Earth¬ 
quake 

H-19 


B.  MULTICHANNEL  FILTERING  OF  VERTICAL  ARRAY  DATA  FROM 
GRAPEVINE  AND  UBO 

During  1965,  deep -well  digital  recording  programs  have 
been  carried  out  at  Grapevine,  Texas,  and  the  Uinta  Basin  Observatory 
(UBO),  It  is  hoped  that  these  data  will  prove  to  be  much  more  useful  for 
vertical  array  evaluation  than  were  any  data  reported  to  data,  since  for  the 
first  time,  3ix  deep  seismometers  were  operated  simultaneously  in  the 
wells.  Instruments  employed  at  Grapevine  were  a  3 -component  seismometer 
at  the  surface  and  six  deep  seismometers  at  3500,  4500,  5500,  6500,  7500, 
and  8500  ft.  Depths  of  the  six  downhole  seismometers  at  UBO  were  3900, 
4900,  5900,  6900,  7900,  and  8900  ft.  In  addition,  the  following  were 
recorded  at  UBO:  1  vertical-component  surface  instrument,  10  vertical- 
component  shallow-buried  seismometers,  2  horizontal -component  surface 
instruments,  outputs  of  10  delay-line  multichannel  filters  operating 
on-line  at  UBO, 


Samples  edited  from  the  Grapevine  recordings  have  been 
described  in  Semiannual  Technical  Report  3  under  Contract  AF  33(657)  - 
12747.  Eighty-seven  samples  were  edited  from  the  first  available  set  of 
UBO  records  including  23  teleseisms,  10  local  events  and  48  noise  samples. 
From  this  set,  the  15  records  listed  in  Table  II - 1  were  selected  for  use 
in  filter  design  and  evaluation. 

At  this  time,  approximately  six  MCF  systems  have  been 
designed  from  and  applied  to  each  set  of  data.  Results  have  not  been 
examined  yet  in  detail  and  will  not  be  discussed  in  this  report,  except 
for  the  following  observation.  Preliminary  study  of  some  UBO  results 
appears  to  confirm  previously  reported  conclusions  that  the  seismic 
noise  at  UBO  is  dominantly  P-waves  between  0.5  and  1.0  cps.  Esti¬ 
mated  signal -to -noise  improvements  for  two  of  the  best  (to  date) 
experimental  UBO  MCF  systems  are  about  5  db  below  0.66  cps,  less  than 
2  db  between  0,66  and  1.66  cps  and  almost  10  db  above  1.66  cps.  These 
figures  probably  are  not  representative  of  the  best  performance  which 
could  be  obtained  with  proper  MCF  design,  but  seem  to  be  reliable 
indicators  of  the  frequencies  at  which  effective  signal  enhancement  can  be 
achieved  at  UBO. 

C.  DECOMPOSITION  OF  UBO  AMBIENT -NOISE  FIELD  INTO 
CONSTITUENT  MODES  * 

From  September  through  October  1964,  3 -channel  vertical 
array  recordings  were  made  at  UBO.  Data  collection  was  described  in 
Semiannual  Technical  Report  2  and  results  of  filter  design  and  evaluation 


^Results  described  in  this  paragraph  have  been  presented  in  'Vertical  Array 
Experiments  at  Uinta  Basin  Seismological  Observatory",  paper  presented  at 
35th  annual  meeting  Society  of  Exploration  Geophysicists,  Dallas,  Texas, 
Nov.  24,  1?65. 
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we"e  reported  in  Semiannual  Technical  Report  3,  both  under  Contract 
AF  33(657) -12747.  Since  observed  MCF  performance  was  much  worse 
than  had  been  predicted  previously  on  a  theoretical  baaip,  *  it  was  deemed 
worthwhile  to  investigate  the  nature  of  the  ambient  noise  at  UBO. 

V 

3y  averaging  values  computed  from  67  recordings,  average 
sets  of  signal  and  noise  autocorrelation  and  crosscorrelation  functions  were 
synthesized  for  a  hypothetical  6 -element  vertical  array  with  receivers  at 
depths  of  200,  900,  2900,  4900,  6900,  and  3900  ft.  For  the  same  array, 
theoretical  sets  of  correlation  functions  have  been  computed  for  P -waves 
t vertically  incident)  and  surface  waves.  These  '?e  displayed  in 
Figure  11-15,  ^ach  set  being  arranged  in  the  sequence  c^,  ep^,  ,  ,  .  , 

<p  6,  CP22,  ....  cp26,  cp33,  .  .  .  cp^  where  the  subscript**  refer  to 

seismometers  in  the  order  of  increasing  depth.  **  The  theoretical  models 
represented  in  Figure  II - 1 5  are  the  following: 

(1)  Vertically  incident  P -waves  with  white  spectrum  in  the  range  of 
0  to  2  cps. 

(2)  Fundamental -mode  surface  waves  with  white  spectrum  in  the 

range  of  0  to  2  cps  •, 

(3)  First  higher-mode  surface  waves  with  white  spectrum  in  the 
range  of  0,  3  to  2  cps 

(4)  Second  higher -mode  surface  waves  with  white  spectrum  in  the 
range  of  0.  5  to  2  cps 

(5)  Uncorrelated  white  noise  with  equal  power  on  all  channels 

Since  the  theoretical  models  were  band-limited  to  the  range 
of  0  to  2.  0  cps,  it  was  decided  to  limit  the  frequencies  in  the  experimental 
correlation  sets  to  the  same  range.  The  upper  traces  in  Figures  11-16 
and  11-17  are  the  band-limited  experimental  correlation  sets  for  noise  and 
signal.  The  aim  of  this  study  is  to  approximate  the  experimental  noise 
correlation  set  by  a  summation  of  a  group  of  correlation  sets  of  known 
or. gin.  The  relative  contributions  ox  the  various  correlation  sets  to  the 
best  summation  may  be  interpreted  then  as  estimates  of  the  relative  con¬ 
tributions  of  the  corresponding  types  of  noise  to  the  observed  noise  field. 


*  Texas  Instruments,  1965,  Array  Research;  horizontal  and  vertical 
arrays  for  teleseismic  signal  enhancement;  UBO  model  theoretical  results: 
Spec.  Rpt.  No.  6,  AFTAC,  Contract  AF  33(657) -12747,  15  Feb. 

##  For  purposes  of  computational  economy,  each  set  of  correlation 
functions  has  been  arranged  in  a  form  analogous  to  a  seismic  trace  and,  in 
this  report,  are  displayed  in  that  form. 
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Figure  II--15.  Theoretical  Correlation  Sets  for  F  -Waves,  Surface  Waves,  and  Incoherent  Noise 


Th*»  simplest  means  of  analyzing  the  experimental  noise 
field  would  be  to  mpute  the  weighted  sum  of  the  theoretical  correlation 
set9  which  provides  the  best  fit  to  the  experimental  correlation  set  on  a 
least -mean-square-error  basis.  This  is  somewhat  undesirable  since  the 
estimated  components  would  be.  constrained  (unreasonably)  to  have  the  same 
spectra  that  were  used  in  computing  their  correlation  functions.  In  order 
to  provide  more  freedom  in  the  estimates,  computations  were  run  to  obtain 
an  optimum  set  of  filters  (one  filter  for  each  theoretic  "I  correlation  set' 
which  couV  be  apj  .ied  to  the  theoretical  correlations. 

The  filters  were  optimized  on  the  basis  that  a  summation  of 
the  filtered  correlation  sets  must  be  the  best  possible  estimate  of  the 
experimental  correlation  set.  Since  the  input  correlations  were  derived 
initially  f'om  an  assumption  of  white  spectra,  and  since  in  the  filtering 
operaf  >n  their  Fourier  t-ansforms  are  multiplied  by  the  amp'itude  responses 
of  the  appropriate  filters,  then  the  power  spectra  associated  with  the 
filtered  correlation  sets  must  be  identical  in  form  to  the  filter  amplitude 
responses.  These  responses  are  interpreted  as  the  best  estimates  of  the 
power  spectra  of  the  contributions  of  the  corresponding  constituents  to  the 
noise  field, 

A  length  of  El  points  (1.44  sec)  was  chosen  for  the  filters. 

The  corres*  onding  frequency  resolution  is  0,695  cps.  If  the  filters  were 
of  infinite  length,  it  would  be  possible  to  make  independent  fits  to  the 
experimental  data  at  every  frequency.  Thi<-,  would  be  almost  equivalent  to 
the  method  used  by  A.  J.  Seriff  in  the  analysis  of  deep-well  data  from 
Juno,  Texas*  (a  "cuiet"  site  similar  to  UBO  in  many  aspects).  In  that 
case,  the  decomposition  was  performed  in  the  frequency  domain.  At 
each  of  four  frequencies,  optimum  weJghting  functions  were  determined  to 
estimate  observed  power  ratios  and  coherences  of  seismic  noise  from 
theoretical  power  ratios  and  coherences  computed  for  P -waves  and  surface 
waves.  These  weights  were  interpreted  as  estimates  of  the  power  densities 
of  the  modes  at  the  corresponding  frequencies.  When  short  optimum 
filters  are  designed  in  the  time  domain,  the  Eolations  are  forced  to  be  smooth 
functions  of  frequency.  Statistical  significance  is  gained  at  the  expense  of 
frequency  resolution. 

The  above  discussion  has  been  concerned  -;ith  the  approximation 
of  the  experimental  noise  correlation  set  by  the  group  of  five  theoretical 
correlation  gets  (Noise  Estimate  1),  The  best  estimate  obtained  in  this 
way  la  shovn  by  the  second  trace  in  Figure  11-16.  The  third  trace  contains 
the  error  set.  — the  difference  between  the  correlations  being  decomposed 
aud  the  best  estimate.  In  addition,  three  other  decompositions  were 
performed: 

♦Shell  Development  Company,  1965,  Seismic  wise  m  deep  boreholes: 

Semiannual  Tech.  Rj  t.  No.  3,  Contract  AF  19(62f5)-2785. 


•  Noi3e  Estimate  2  — as  above,  but  with  the  theoretical 
P -wave  correlation  set  replaced  by  the  experimental 
signal  set 

•  Signal  Estimate  1  — approximation  of  the  experimental 
signal  correlation  set  by  a  combination  of  the  five 
theoretical  sets 

•  Signal  Estimate  2  —approximation  of  the  experimental 
signal  correlation  set  by  a  combination  of  the  theoretical 
sets  for  P-waves  and  'ncoherent  noise  and  the  experimental 
noise  correlation  set 

The  optimum  summation  traces  and  the  associated  error 
traces  are  shown  in  Figures  11-16  and  11-17.  The  estimated  power  spectra 
(iae.,  the  optimum  filter  responses)  are  plotted  in  Figure  71-18. 

The  results  indicate  that  P-waves  f^rm  a  major  part  of  the 
UBC  ambient -noise  field  below  1  cps  and  that  the  fundamental  Rayleigh  mode 
is  much  more  important  at  the  surface  than  any  of  the  higher-order  modes. 
This  finding  is  in  general  agreement  with  the  results  obtained  for  Juno, 

Texas,  from  which  it  was  estimated  that  P-waves  and  surface  waves  were 
of  approximately  equal  strength  at  0.25  and  0.  50  cps  and  thac  P-wave  power 
was  about  -5  db  relative  to  surface-wave  power  at  1  cps.  Further  support 
for  the  hypothesis  of  the  existence  of  strong  P-wave  noise  at  Juno  was 
provided  by  the  observation  of  noise  bursts  with  velocities  greater  than 
8  km/sec.  The  UBO  1-cps  wavenumber  spectra  (Semiannual  Technical 
Report  3,  Contract  AF  33(657) -12747  indicate  that  about  20  perr  nt  of  the 
ambient  noise  is  P-waves  and  that  the  only  other  important  contributor 
is  the  fundamental  mode.  Unfortunately,  the  UBO  array  is  not  sufficiently 
large  to  permit  useful  resolution  in  wavenumber  space  at  crequenoies  vary 
much  lower  than  1  cps. 

In  both  signal  and  noise  predictions,  'It  was  found  that  the 
residual  error  was  less  for  Estimate  2  (obtained  xrom  a  mixture  of 
experimental  and  theoretical  correlation  rets)  than  for  Estimate  1 
(obtained  from  theoretical  sets  alone).  Thi3  difference  is  most  pronounced 
at  low  frequencies,  especially  below  0.  25  cps;  thus,  the  low-frequency  parts 
of  the  two  experimental  set,  resemble  each  other  much  more  closely  than 
any  of  the  theoretical  models.  Since  low-frequency  3/N  is  very  poor  in  the 
signals  recorded  at  UBC.  it  is  reasonable  that  "signal"  statistics  and 
"noise"  statistics  should  appear  similar.  Lack  of  agreement  with  theoretical 
values  probably  is  uue  to  the  failure  to  correct  for  seismometer  responses 
at  frequencies  other  than  1  cps.  Whet  the  process  was  forced  *o  choose 
among  theoretical  monels,  the  results  bad  cated  high  u-w.,ve’  to -surface - 
wave  power  ratios  at  low  frequencies  in  both  cases.  In  spite  of  disagreements 
between  theoretical  and  experimental  statistics,  the  process  was  able  to 
make  ».  clear  choice  in  favor  of  P-wave  :  .ns a. 


D„  INVESTIGATIONS  OF  AMBIENT  P-WAVE  NOISE 


Results  obtained  from  a  number  of  recording  sites  indicate 
the  existence  of  a  P-wave  constituent  in  the  ambient  seismic  noise.  It 
ha3  been  suggested  that  the  P-wave  noise  is  caused  by  a  steady  background 
of  P -waves  emerging  from  the  mantle  with  uniform  power  on  a  worldwide 
scale.  It  would  follow  that  the  only  spatial  variations  observed  at  the 
surface  should  be  thos  »  which  are  due  to  differences  in  crustal  transmission 
characteristics.  Such  differences  would  not  be  expected  to  produce 
variations  in  observed  power  of  more  than  a  factor  of  two  {assuming  good 
coupling  of  sensors  to  solid  unweathered  rock). 

Multichannel  filter  systems  employed  in  conjunction  with  3 -km 
(approximately)  arrays  at  "noisy"  sites,  such  as  Wichita  Mountains  and 
Cumberland  Plateau,  have  been  found  to  be  very  effective  for  wide -band  ambient  - 
noise  rejection.  At  "quiet"  sites,  such  as  Uinta  Basin  and  Tonto  Forest, 
significant  signal -to -noise  improvement  has  not  been  obtainedin  the 
frequency  range  0.  5-1.2  cps.  These  findings  are  consistent  wdth  the 
hypothesis  of  a  constant  level  of  P-wave  noise,  since  it  appears  that 
0.5-1.  2  cps  noise  outputs  of  processors  applied  to  small  (less  than  5  km) 
arrays  are  restricted  by  a  lower  limit  which  is  much  more  uniform  than 
the  input -noise  power  levels.  If  information  below  1.2  cps  is  essential, 
then  the  only  knov/n  method  of  penetrating  this  barrier  is  the  use  of  large 
arrays  such  as  the  200 -kro  LASA  array  in  Montana.  With  such  an  array, 
it  should  be  possible  to  use  processing  systems  tuned  to  pass  only  P -waves 
originating  from  specified  regions  of  the  globe.  It  is  clear  that  P-wave  noise 
must  be  a  considerat  n  of  prime  importance  in  the  design  and  location 
choice  of  future  array  istallations. 

Figure  11-19  shows  ambient -noise  power  spect  a  obtained 
from  five  recording  sites.  It  is  believed  that  these  curves  should  be 
related  to  P-wave  noise  level  as  follows: 

0  UBO:  This  spectrum  was  obtained  by  applying  the 
computed  fractional  contribution  of  P -waves  to  the 
total  noise  powur  spectrum  for  a  single  noise  sample 
and  smoothing.  This  appears  to  be  a  reasonable  estimate 
of  P-wave  noise  power  in  the  neighborhood  of  1  cps,  but 
probably  is  too  high  at  lower  frequencies  because  of 
the  3trong  microseismic  peak  in  the  original  data  and 
the  forced  smoothness  of  the  estimated  P-wave  fraction. 

•  TFO:  This  curve  is  the  spectrum  of  a  noise  sample  as 
recorded  by  a  eingle  seismometer.  Since  the  0.5-1. 2  cps 
noise  at  TFO  is  assumed  to  be  almost  entirely  P-waves, 
the  plotted  values  should  be  only  slightly  higher  than 
the  P-wave  spectrum  in  this  range. 


RELATIVE  TO  l(niii)  /CF&  AT  1  CPS 


Figure  11-19.  Ambient  Noise  Spectra  from  Five  Array  Site3 
(Not  Corrected  for  Frequency  Response  of 
J -M  Seismometers) 


•  CPO:  This  spectrum  was  computed  for  a  noise  sample 
after  it  was  nassed  through  a  signal  enhancement  filter 
(isotropic  velocity  filter).  It  is  believed  that  the 
output  noise  is  due  mainly  to  P -waves  below 

1.0  cps  and  surface  waves  or  uncorrelated  noise 
above  1. 0  cps. 

•  PI2WY:  The  curve  represents  the  lower  envelope  of 
spectra  computed  for  the  six  PI2WY  records  (paragraph 
A).  If  a  constant  P-wave  noise" level  is  assumed,  then 
the  total  noise  power  always  must  be  at  least  as  great 
as  that  constant  level.  Therefore,  the  curve  plotted 
must  be  an  upper  limit  to  the  P-wave  spectrum. 

The  poor  performance  of  MCF  systems  designed  for 
this  site  implies  that  the  P-wave  spectrum  is  not  more 
than  3  db  below  the  limit  in  ie  0.  5-1.2  cps  range. 

•  WMO:  A  simple  summation  of  the  outputs  of  the  10- 
seismometer  array  was  obtained  and  a  power  spectrum 
was  computed  for  a  single  noise  sample.  Surface -wave 
rejection  was  less  effective  here  than  in  the  case  of  the  CPO 
processor  and  so  this  curve  can  only  be  interpreted  as  an 
upper  bound  on  the  possible  P-wave  noise  spectrum. 

If  the  low-frequency  part  of  the  UBO  curve  is  discarded,  then 
the  agreement  becomes  very  good.  It  should  be  noted  that  the  P-wave 
noise  power  density  reported  by  Seriff  for  Juno,  Texas,  is  approximately 
-13  db  relative  to  1  mM^/cps  at  1  cps.  This  independent  estimate  is  also 
in  good  agreement.  The  curves  of  Figure  II - 1 9  are  not  corrected  for  the 
J-M  seismometc ,  response.  Application  of  this  correction  to  an  average 
P-wave  spectrum  deduced  from  those  curves  yields  the  absolute  P-wave 
power  spectrum  in  Figure  11-20. 

Experimental  evidence  supports  the  hypothesis  of  P-wave 
noise  with  a  uniform  (within  a  factor  of  2)  worldwide  power  level.  An 
estimated  figure  of  0.05  mn2/CpS  at  1  cps  appears  to  be  accurate  to  better 
than  an  order  of  magnitude.  The  only  known  method  for  enhancing  desired 
P-wave  signals  in  the  presence  of  such  a  background  of  unwanted  P-wave 
noise  is  to  improve  the  wavenumber  resolution  of  an  array  by  increasing 
its  size  to  hundreds  of  kilometers,  if  signal  and  noise  both  occupy  the 
same  frequency  band.  The  shape  of  the  estimated  spectrum  of  P-wave 
noise  implies  that  some  advantage  might  be  gained  by  a  concentration  of 
effort  in  the  neighborhood  of  2  cps  instead  of  1  cps.  It  is  felt  that 
vertical  arrays  and  small -aperture  horizontal  arrays  may  well  offer  useful 
advantages  at  frequencies  above  1.  5  cps.  Such  advantages  could  probably  be 
exploited  more  effectively  by  vertical  arrays  with  smaller  receiver  spacings 
than  have  been  used  to  date  . 
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Figure  11-20.  Estimate  o£  Absolute  P-Wave  Noise  Power  Spectrum 
(After  Instrument  Response  Correction) 


The  spectrum  presented  in  Figure  II -2C  ig  admittedly  a  very 
crude  estimate.  It  is  planned  to  attempt  to  refine  the  mode  separation 
results  for  UBO  and  also  to  apply  similar  techniques  to  data  from  other 
vertical  arrays.  It  should  be  possible  to  gain  additional  useful  information 
from  a  study  of  ambient -noi3e  wavenumber  spectra  which  have  been  computed 
for  a  number  of  horizontal  array  installations. 
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SECTION  III 


STUDY  OF  TELESEISMS  RECORDED  AT  CFO 
A.  INTRODUCTION 

Three  ensembles  of  teleseisms  recorded  by  multiple  array 
processor  (MAP)  systems  at  Cumberland  Plateau  Observatory  (CPO) 
have  been  and  are  being  studied.  Almost  all  attention  has  been  restricted 
to  the  records  output  by  isotropic  processors  which  weru  designed  to  reject 
low-velocity  noise  and  preserve  high-velocity  signals  with  a  minimal 
dependence  on  azimuth.  These  records  contain  relatively  clean  body-wave 
signals  since  both  ambient  noise  and  signal -generated  noise  (energy  scattered 
into  surface  waves  in  the  vicinity  of  the  receiver)  are  attenuated  by  the 
processors. 


Figure  III  - 1  shows  an  example  of  the  processor  performance. 
For  a  representative  event,  pov/er  spectra  were  computed  for  an  interval 
containing  the  signal  and  an  interval  p^eceeding  the  signal  onset.  Between 
0.  75  and  2  cps,  a  signal -to-noise  improvement  of  about  10  to  15  db  is  observed 
for  the  processor  relative  to  a  single  seismometer.  Since  the  estimated 
signal  spectra  are  actually  spectra  of  signal  plus  noise,  the  improvement 
is,  in  fact,  somewhat  greater  than  it  appears. 

Objectives  of  the  study  are  described  in  the  following 
paragraphs.  In  some  cases,  only  brief  outlines  of  the  results  will  be 
given,  since  more  detailed  accounts  may  be  found  in  Semiannual  Technical 
Reports  1,  2  and  3  under  Contract  AF  33(657) -1 2747. 

B.  CRUSTAL  REVERBERATION  DECONVOLUTION 

An  impulsive  P -wave  signal  emerging  from  the  mantle 
encounters  a  series  of  reflecting  interfaces.  The  resulting  reverberations 
of  the  signal  have  the  effect  of  a  frequency  filter  which  converts  the  impulse 
into  a  wavelet  of  theoretically  infinite  duration.  In  practice,  the  time  span 
of  the  wavelet  is  finite  since  the  amplitude  cannot  remain  above  detection 
threshold.  It  would  be  desirable  to  design  an  inverse  filter  which  would 
operate  on  the  observed  signal  in  such  a  way  as  to  remove  crustal  effects  by 
transforming  the  waveform  back  to  an  impulse. 

As  reported  in  the  Semiannual  Technical  Reports  2  and  3, 
attempts  have  beer:  made  to  apply  dereverboration  processing  to  the  events 
of  all  three  ensembles.  Filter  design  was  based  upon  signal  auto¬ 
correlations  computed  from  data  samples  containing  the  P -phase. 
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These  correlations  depend  on  crustal  effects  to  a  much  lesser  degree  than 
on  incident  signal  spectra,  which  depend  on  both  source  spectra  and  mantle 
transmission  effects  (including  anelastic  absorption).  Thus,  the  major  effect 
cf  a  deconvolution  filter  is  to  compensate  for  the  incident  signal  spectrum 
which  is  assumed  to  be  severely  band-limited,  but  not  to  exhibit  fine 
structure.  However,  if  the  filter  is  sufficiently  long,  the  fine  structure 
reduced  by  crustal  effects  may  also  be  compensated  for, 

Unless  the  span  of  the  deconvolution  filter  is  more  than  a  few 
seconds,  differences  in  emergence  angle  will  not  be  important.  For  example, 
Figure  III -2  shows  that  some  fine  structure  (probably  associated  with  Moho 
reflections)  persis  ted  even  when  power  spectra  for  60  teleseisms  of  world¬ 
wide  origin  were  averaged.  Thus,  it  is  possible  to  design  an  "average" 
filter  which  will  deconvolve  an  "average"  teleseismic  signal.  It  was  shown 
in  the  case  of  Ensemble  I  that  an  ensemble  average  deconvolution  filter 
was  almost  as  effective  for  waveform  contraction  as  filters  designed 
individually  for  each  event. 

Although  some  crustal  dereverbera.tion  may  be  achieved  with 
deconvolution  filters,  it  must  be  borne  in  mind  that  the  filter  response 
cannot  be  interpreted  directly  as  the  inverse  of  the  *rustal  response  since 
the  most  important  effect  is  compensation  for  the  incident  signal  spectrum. 
Since  signal -to-noise  ratio  is  generally  poor  wherever  signal  power  is  low, 
this  compensation  will  tend,  in  general,  to  amplify  noise  more  than  signal. 

It  is  therefore  advisable  to  apply  band -limiting  filters  to  the  deconvolved 
data  in  order  to  reject  noise  in  the  frequency  ranges  where  signal  power  is 
low.  If  these  filters  have  smooth  frequency  responses,  compensation  for 
fine  structure  will  have  been  achieved. 

C.  DEPTH  OF  FOCUS  DETERMINATION 

Focal  depths  of  earthquakes  are  determined  most  often  by 
least-mean -square  fics  to  observed  arrival  times  of  P.  An  alternative 
method  depends  upon  the  difference  in  observed  arrival  times  of  P  and  pP. 

The  pP  -  P  method  probably  is  more  accurate,  but  requires  that  pP  be 
observed  and  identified  with  certainty.  The  value  of  multichannel  velocity 
filtering  of  seismometer  array  outputs  has  been  demonstrated  in  that  it  has 
been  possible  sometimes  to  identify  pP  in  the  output  trace  of  the  CPO 
processor  where  it  was  not  observable  in  single -seismometer  records. 
However,  in  other  cases,  pP  could  not  be  seen  even  when  its  predicted 
arrival  time  was  favorable  for  observation  (i.  e. ,  when  USC&GS  estimated 
focal  depth  indicated  pP  should  arrive  after  P  has  died  out,  but  before  end 
of  the  record). 
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Average  Power  Spectrum  of  60  Ensemble  II  Teleseisms 


It  has  been  suggested  that  it  might  be  possible  to  improve  pP 
detectibility  by  deconvolving  the  processor  output  record  or  by  computing 
autocorrelations  for  either  the  raw  t::ace  or  the  deconvolved  trace.  None 
of  these  techniques  was  found  to  offer  improvements.  This  failure  has 
been  attributed  to  differences  between  the  waveiorms  of  P  and  pP  from  the 
same  events. 

It  must  be  concluded  that  pP  -  P  analysis  is  of  limited  value 
in  focal  depth  estimation  for  earthquake /explosion  discrimination.  A  well- 
documented  pP  phase  offers  excellent  proof  that  focal  depth  was  greater 
than  a  few  kilometers.  However,  several  deep -fo '.us  (USC&GS  estimated 
depth  >100  km)  events  did  not  produce  observable  pP  arrivals,  A  possible 
explanation  is  asymmetry  of  the  radiation  pattern  at  the  source.  Regardless 
of  the  reasons  for  pP  absence,  it  must  be  concluded  that  such  absence 
cannot  be  accepted  as  evidence  of  near-surface  origin. 

D.  STUDY  OF  CPO  CRUSTAL  STRUCTURE 

Two  important  effects  must  be  considered  when  deduction  of 
crustal  structure  beneath  CPC  from  analysis  of  teleseismic  records  is 
attempted.  Since  crustal  response  depends  on  signal  angle  of  emergence, 
it  is  not  meaningful  to  estimate  a  single  average  response  function  for  the 
complete  set  of  teleseisms.  In  the  case  of  Ensemble  I,  ?3  events  from  a 
small  range  of  distances  were  available.  By  forming  seven  different  sub¬ 
groups  and  computing  average  autocorrelations  (Semiannual  Technical 
Report  3),  it  was  possible  to  find  "events”  in  the  autocorrelation  functions 
which  could  be  interpreted  as  evidence  of  major  reflecting  horizons  below 
CPO.  In  Figure  III -3,  it  can  be  seen  that  such  "events"  no  longer  appear 
when  ensembles  of  worldwide  origin  are  averaged.  The  one  coincidence 
indicated  at  18-sec  lag  is  probably  not  physically  significant. 

Also,  observed  spectra  can  be  interpreted  as  crustal  response 
functions  only  after  they  have  been  corrected  for  spectral  content  of  the 
original  signal  incident  at  the  base  of  the  crust.  Figure  III  — i  shows  average 
autocorrelations  of  all  Ensemble  II  and  III  events  between  60°  and  75°. 
Evidence  of  reflections  has  been  smeared  cut  beyond  recognition,  probably 
because  of  differences  in  signal  spectra.  Figure  III  - 5  illustrates  the  great 
diversity  observed  in  teleseisms  from  a  small  range  of  distances.  Power 
spectra  have  been  computed  for  18  Ensemble  II  events  between  ?0°  and 
30°.  It  can  easily  be  demonstrated  that  these  differences  are  associated  with 
the  signals,  since  noise  power  spectra  (Figure  III -6)  computed  for  intervals 
just  preceding  each  event  show  excellent  time -stationarity.  It  is  believed 
that  observed  differences  in  signals  from  nearly  the  same  distance  may  be 
explained  by  near -source  effects.  It  should  be  possible  to  average  out 
these  effects  to  some  extent  by  the  use  of  large  numbers  of  teleseisms. 

Text  cont'd  page  III- 10 
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In  Ensembles  II  and  III,  it  is  not  possible  to  obtain  such  a 
concentration  of  foci.  It  appears  that  the  most  effective  method  o£ 
eliminating  source  effects  might  be  to  discard  all  shallow-iocus  events. 

As  an  example,  power  spectra  for  seven  deep-focus  events  {11-21  and  III-23 
through  111-30)  are  shown  in  Figure  III -7;  it  has  not  been  possible  to  find 
such  striking  resemblances  in  similar  groups  of  shallow -focus  events.  In 
Figure  III -8,  autocorrelations  of  four  teleseisms  (III -27  through  III  -  30) 
show  four  "events"  which  are  reproducible  throughout  the  set.  A  the  rough 
study  of  the  available  data  has  failed  unfortunately  to  turn  up  any  ether  sets 
which  display  such  similarities .  The  four  earthquakes  which  contributed  to 
Figure  III -8  occurred  at  distances  of  30.9°,  30.  f  3  ,  31,33  ,  and  3l.0l 
with  azimuths  between  155.1°  and  155.6°  and  focal  depths  between  140  and 
176  km.  This  set  was  the  only  available  set  of  deep -focus  earthquakes 
from  such  a  small  area. 

A  major  difficulty  in  interpreting  the  autocorrelations  and 
power  spectra  in  terms  of  crustal  structure  arises  from  the  limited  band¬ 
width  of  the  signals  being  studied.  Since  most  signals  have  useful  band- 
widths  of  about  1  cps,  the  finest  resolution  which  can  be  obtained  is  about 
1  sec  of  time  or  a  few  kilometers  of  depth.  A  suggested  model  for  the 
structure  under  CPO  is  reproduced  in  f  igure  III-9.*  Since  the  two-way 
travel  time  through  the  entire  sedimentary  section  is  only  1,05  sec  for 
this  model,  it  is  impossible  to  learn  anything  about  detailed  structure  in 
the  sediments  and  very  difficult  to  determine  even  their  total  thickness. 

Ideally,  autocorrelation  functions  should  exhibit  peaks  at  time  lag 3 
corresponding  to  the  differences  in  arrival  times  of  direct  and  reflected 
signals.  The  correlations  in  Figures  III  -  3  and  III -8  are  computed  from 
signals  with  phase  velocities  of  2?.  5  km/sec  (Kurile  Islands,  85  )  and 
12.6  km/sec  (Colombia  Venezuela,  31°).  The  difference  in  phase  velocities 
should  be  reflected  by  shifts  of  the  peaks  in  the  autocorrelation.* . 

Some  events  that  are  multiply  reflected  and  are  most  likely 
to  contain  significant  energy  are  illustrated  in  Figure  III- 10.  The  time  lags 
of  these  events  relative  to  the  direct  P-wave  are  given  in  Table  III  - 1  for 
the  crustal  model  of  Figure  III  —  9 .  Some  features  in  the  autocorrelations 
have  been  sel  ‘cted  as  evidence  of  reflected  arrivals,  but  it  is  difficult  to 
decide  upon  the  proper  lag  times  which  should  be  assigned  to  such  features. 
Ideally,  a  reflected  arrival  should  give  rise  to  a  reproduction  of  the  center 
part  of  the  autocorrelation,  symmetrical  and  centered  at  the  lag  time  of  the 
arrival.  Therefore,  it  is  preferable  to  seek  center  times  rather  than  onset 
times  of  the  "events."  Since  most  >f  the  "events"  selected  are  unsymmetrical, 
center  times  are  difficult  to  estimate.  The  added  complication  of  limited 

^Foreman,  J. ,  ?  964,  Some  theoretical  calculations  on  a  model  of  the  crust 
under  the  Cumberland  Plateau  Observatory:  Texas  Instruments  (Unpublished). 

Text  cont'd  page  111-15 
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Figure  III ~ 9 -  Model  of  Crustal  StructureNear  McMinnille,  Tennessee 


Figure  III-10.  Travel-Paths  of  Reflected  Arrival 


signal  bandwidth  implies  that,  opt  mum  resolution  is  between  1  and  2  sec. 


Table  1II-1 

COMPUTED  LAG  TIMES  OF  REFLECTED  ARRIVALS 


—1 

Signal 
Phase 
Velocity 
fkm/  sec) 

Arrival  A 

Arrival  B 

Lag  Time  (sec) 

Arrival  E 

Arrival  F 

Arrival  C 

Arrival  D 

12.6 

10.  0 

6.  0 

**•  o 

9.0 

5.0 

14.  0 

22.  5 

11.0 

6.  5 

4.  5 

10.9 

5.  5 

15.  5 

03 

11.5 

6.  8 

4.  7 

10.4 

5.7 

16.2 

Estimated  times  for  the  first  two  "events'1  in  Figure  III “3 
are  6.0  sec  and  10.5  sec.  These  probably  represent  events  B  and  A  with 
22.  5-km/sec  phase  velocity.  If  this  identification  is  correct  and  the  lag 
times  have  been  estimated  accurately,  then  the  depths  of  the  Conrad  and 
Mohorovicic  discontinuities  as  given  in  Figure  III-  )  are  too  large  by 
about  2  km.  However,  the  dominant  near -surface  reflector  could  be  an 
interface  within  the  sedimentary  section  or  the  base  of  the  sediments,  in 
which  case  the  predicted  Conrad  and  Mohorovicic  depths  might  even  be  too 
small  by  about  2  km.  The  "event'1  founa  near  15  sec  could  be  any  one  of  a 
large  set  of  multiply -reflected  arrivals,  but  event  F  seems  to  be  the  most 
likely  arrival.  The  ambiguity  in  interpretation  remains,  since  the  apparent 
error  in  lag  time  can  be  removed  either  by  decreasing  the  Conrad  a'nd 
Mohorovicic  depths  or  by  assuming  the  dominant  near -surface  reflector  to 
be  the  base  of  the  Devonian. 

In  Figure  III >8,  "even's  1  have  been  found  at  4.  2  and  12.  ? 
sec  for  a  phase  velocity  of  12.6  km/sec.  The  major  peaks  In  Figure  J.II-7 
occur  at  intervals  of  0,28  cps,  implying  a  reflecting  system  with  a  lag  time 
of  3.6  sec.  The  differences  in  the  lag  times  found  for  22.  6  k/n/see  and 
12.6  km/sec  are  much  greater  than  would  be  expected,  unle  ss  different 
reflecting  interfaces  predominate  in  the  two  cases.  It  is  difficult  to  conceive 
of  how  a  model,  consistent  with  Figure  III  —  9  ,  could  give  rise  to  only  the 
"events"  seen  in  P'igure  III-  3.  Since  the  four  earthquakes  contributing  to 
Figure  III- 8  were  so  similar,  it  may  be  that  the  observed  coincidences  can 
be  attributed  to  source  effects. 
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In  summary,  the  result  j  suggest  that  the  assumed  crustal 
model  for  CPO  is  not  '.n  er  rnv  by  more  than  1  or  2  km,  but  the  results  are 
not  of  sufficient  quality  to  permit  more  definite  conclusions  at  this  time. 

3.  INVESTIGATION  OF  PROPAGATION  MECHANISMS 

An  attempt  is  being  made  to  use  the  high-quality  teleseism 
ipfords  '  orn  the-  CPO  processor  do  examine  such  subjects  as  inelastic 
absorption  in  the  earth's  interior  and  the  dependence  of  crustal  response  on 
3ourr.‘!  range  and  azimuth.  Therefore,  moi  e  than  100  different  groups  of 
e /ents  have  been  chosen  from  Ensembles  II  an:'.  Ill,  and  autocorrelation 
functions  ana  powe’-  spectra  have  been  computed  for  each  subensemb’e. 

These  data  have  not  been  studied  yet  in  great  detail.  A 
preliminary  examination  suggests  that  because  of  the  large  number  of 
variables  involved  (e.g.  ,  focal  depth,  geographical  location  of  source, 
source-to-receiver  distance  and  azirmith,  geological  source  environment, 
source  spectrum,  magnitude,  continental  or  oceanic  nature  of  early  part 
of  travel  path,  source  radiation  pattern),  the  number  of  available  recordings 
is  too  small  to  provide  useful  statistical  significance. 
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SECTION  IV 


A  SUMMARY  OF  CPO-MCF  SYSTEMS 

A.  INTRODUCTION  AND  CONCLUSIONS 

This  is  a  summary  o £  work  which  soon  will  be  published  as 
a  special  report  under  the  title  "A  Study  of  Multichannel  Filter  Systems 
Incorporating  Statistical  Gain  Fluctuation  Designed  for  Cumberland  Plateau 
Observatory." 


MCF  systems  which  were  indicated  to  have  had  gain  inequality 
problems  were  redesigned,  incorporating  statistical  gain  fluctuation.  (Under 
certain  circumstances,  gain  inequalities  can  cause  multichannel  filter 
systems  to  separate  noise  and  signal  on  the  basis  of  these  inequalities 
rather  than  on  true  differences  in  spatial  organization.  When  this  happens, 
a  system  is  said  to  have  gain  inequalization  problems.)  The  counterpart 
systems  then  were  compared  in  detail.  This  study  also  had  a  bearing  on  the 
evaluation  of  the  partial  arrays  at  CPO^  in  that  the  multichannel  filter  systems 
designed  on  the  various  partial  arrays  were  thought  to  have  been  somewhat 
affected  by  gain  inequalization  problems. 

The  results  of  this  study  indicate  that 

»  The  addition  of  statistical  gain  fluctuation  was  a  good 
solution  to  the  gain  inequalization  problem. 

•  The  addition  of  statistical  gain  fluctuation  essentially 
prevents  high  gains  in  the  frequency  response  of  the  filters 
for  individual  channels  which  appear  to  be  necessary  for 
MCF  systems  to  operate  cn  the  basis  of  gain  inequalities. 

#  The  conclusions  about  the  effectiveness  of  the  partial 
arrays  at  CPO  which  previously  were  published  ^ 
appear  to  be  vaJid, 

L.  CPO-i  F  SYSTEMS  WITH  GAIN  INEQUALIZATION  PROBLEMS 

'linen  MCF  systems  are  designed  in  the  time  domain  using 
equalized  signal  correlations  and  measured  noise  correlations,  the  systems 
may  middesign.  When  the  seismometers  are  gained  unequally,  the  MCF 
system  can  use  the  inequalities  in  the  noise  correlations  as  compared  to 
the  equalized  signal  correlations  to  separate  the  noise  and  the  artificial 
signal.  This  has  been  found  to  most  strongly  affect  the  filters  at  low 
frequencies  where  the  velocity  differences  between  signal  and  noise  are 
less  significant  and  the  noise  is  highly  coherent.  The  system  so  designed 
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could  be  far  from  an  optimum  system  for  actual  signal -noise  separation, 

MCF  systems  IP  13,  IP  23  through  IP  26  were  designed  for 
CPO.  These  systems  showed,  in  varying  degrees,  symptoms  characteristic 
of  MCF  systems  which  use  seismometer  gain  differences  to  give  artificial 
signal -to -noise  improvement. 

C°0  IP  13  was  designed  as  a  5 -channel  system  to  extract 
the  signal  at  the  center  seismometer  of  the  array  (Figure  IV-1).  The  5 
channels  were  the  center  seismometer,  the  sum  of  the  3  seismometers  of 
ling  2,  the  sum  of  the  3  seismometers  of  ring  3,  the  sum  of  the  6  seis¬ 
mometers  on  ring  4,  and  the  sum  of  the  6  seismometers  on  ring  5.  IP  13 
was  designed  in  the  time  domain  using  CPO  noise  sample  C  °  and  an 
infinite  velocity  signal  model. 

It  has  been  shown  that  IP  13  was  affected  adversely  by 
inequalization  of  the  seismometers  in  particular,  this  system  rejected 
seismic  noise  as  recorded  through  the  seismometer -PTA  system  by  as 
much  as  23  db  at  0.4  cps,  passed  perfectly  equalized  synthetic  signals  with 
very  little  signal  distortion  and  rejected  actual  signals  which  were  recorded 
through  the  seismometer-PTA  system.  At  low  frequencies,  IP  13  exhibited 
the  characteristics  of  systems  which  are  designing  on  the  basis  of  seis¬ 
mometer  gain  inequalities:  high  gains  in  the  filters  on  individual  channels, 
K-space  power  response  of  the  MCF  system  not  amendable  to  rejection  of 
the  spatial-organized  noise.  (See  Figures  IV-2  and  IV-11,  presented  later.) 

The  MCF  systems  IP  23,  IP  24,  IP  25,  and  IP  26  were 
designed  to  test  the  effectiveness  of  partial  arrays  at  CPO.  These  systems 
have  been  discussed  in  detail  in  a  previous  publication  ■*.  They  were 
designed  to  extract  signal  at  the  center  seismometer  with  IP  23  using  5 
channels,  the  center  seismometer  and  each  summed  ring;  IP  24  using  4 
channels,  the  center  seismometer  and  the  summed  rings  3,  4  and  5; 

IP  25  using  3  channels,  the  center  seismometer  and  the  summed  rings  4 
and  5;  IP  26  using  2  channels,  the  center  seismometer  and  the  summed 
ring  5  (Figure  IV-1).  These  systems  were  designed  in  the  time  domain 
using  an  infinite  velocity  signal  model. 

These  systems  were  thought  to  have  been  affected  by 
seismometer  gain  differences  because  the  systems  showed  differences  in 
their  K-space  responses  at  0.  5  cps.  Also,  the  low-frequency  characteristics 
of  the  individual  channel  filters  implied  seismometer  gain  differences 
(p.  5)  The  symptoms  of  gain  in  equalization  problems  for  these  systems 
were  less  severe  then  those  of  IP  13. 
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C.  DESIGN  OF  COUNTERPART  MCF  SYSTEMS  WITH  GAIN  FLUCTUATION 

Using  the  method  of  incorporating  statistical  gain  fluctuation  u, 
all  of  the  multichannel  filter  systems  mentioned  above  were  redesigned. 

The  effect  of  this  gain  fluctuation  was  to  add  spatially  random  signal  to  the 
multichannel  filter  design  problem.  IP  13  was  redesigned  as  IP  22;  IP  23 
through  IP  26  were  redesigned  as  IP  33  through  IP  36.  IP  22  differed 
from  IP  13  in  that  the  method  of  Incorporating  statistical  gain  fluctuation 
resulted  in  all  signal  autocorrelations  being  scaled  up  and  that  the  time  domain 
filters  were  31  points  long  for  IP  22  as  compared  with  25  points  for  IP  13. 

The  MCF  systems  IP  33  through  IP  36  differed  from  their  counterparts 
IP  23  through  IP  26  only  ir.  that  statistical  gain  fluctuation  was  incorporated 
into  the  redesigned  systems. 

D.  EFFECTS  OF  INCORPORATING  STATISTICAL  GAIN  FLUCTUATION 

The  method  of  adding  statistical  gain  fluctuation,  in  effect, 
adds  some  spatially  uncorrelated  signal.  The  addition  of  spatially 
uncorrelated  signal  makes  high  gains  on  any  channel  unprofitable.  If 
a  large  gain  is  applied  to  one  channel,  it  must  be  compensated  for  on 
another  (to  maintain  signal  response).  The  spatially  uncorrelated  signal 
power  will  be  gained  up  by  both  channels,  with  no  cancellation  occurring, 
which  results  in  a  large  contribution  to  the  mean-square-error  (MSE).  The 
magnitude  of  this  effect  would  depend  upon  signal  power  at  any  frequency. 

High  gains  on  individual  channels  seem  to  be  necessary  ^  for  a  MCF 
system  which  suppresses  noiee  mostly  on  the  basis  of  inequalized  noise 
discriminated  from  equalized  signal.  The  method  of  incorporating  gain 
fluctuation  makes  these  large  gains  unprofitable,  thus,  forcing  the  MCF 
system  to  separate  signal  and  noise  on  the  basis  of  correlation  differences 
which  are  the  result  of  the  spatial  organization  of  the  signal  and  noise. 

Spatially  uncorrelated  signal  powrer  tends  to  force  the  MCF 
system  to  use  the  center  seismometer  for  signal  preservation.  Since  the 
system  wishes  to  output  an  approximation  of  the  signal  as  seen  on  the  center 
seismometer,  it  is  obvious  that  using  any  other  channel  to  estimate  the 
signal  will  increase  the  MSE.  To  produce  a  least  MSE  system  at  each 
frequency,  the  MCF  system  must  weigh  the  increased  signal  distortion, 
which  results  from  using  the  center  seismometer  to  reject  ncise  (and  thus, 
necessarily  other  channels  to  estimate  the  signal),  against  the  decreased 
noise  power  output  possible  by  using  the  center  seismometer  to  reject 
noise.  This  weighing  at  a  particular  frequency  (in  terms  of  cost  in  MSE) 
would  be  effected  by: 


•  The  gains  applied  to  the  ring -summed  channels 


e  The  signal  power 

By  comparing  the  frequency  responses  of  the  individual 
channels  for  the  counterpart  systems  IP  13  with  IP  22,  IP  23  with  IP  33, 

IP  24  with  IP  34,  IP  25  with  IP  35,  IP  26  with  IP  36  (Figures  IV-2  through 
IV-10),  it  is  readily  seen  that  the  systems  with  gain  fluctuation  do  not 
show  the  high  gains  at  low  frequencies  and  are  making  heavier  use  of  the 
center  seismometer  to  pass  the  signal.  The  response  of  channel  I  {center 
seis)  io  smoothed  for  the  systems  with  gain  fluctuation.  These  differences 
are  more  pronounced  as  the  number  of  channels  increases.  The  systems 
using  more  channels  can  better  afford  (in  terms  of  MSE)  to  partially  give, 
up  the  center  seismometer  and  still  suppress  the  noise.  IP  36,  which  has 
only  2  channels,  cannot  possibly  give  up  one  of  them  to  any  large  degree 
and  still  effect  noise  suppression. 

The  notably  greater  dependence  of  the  gain  fluctuation  systems 
on  the  center  seismometer  to  preserve  signal  at  low  frequencies  is  the 
result  of  larger  gains  on  the  ring-summed  channels  and  more  signal  power 
at  these  low  frequencies. 

At  frequencies  where  the  gain  inequalization  problem  strongly 
affected  the  multichannel  filter  systems,  the  K-space  responses  of  tne 
counterpart  eystems  with  gain  fluctuation  are  notably  different.  At 
frequencies  which  were  not  affected  by  gain  equalization  problems,  the 
K-space  responses  of  the  counterpart  systems  generally  are  similar.. 

At  frequencies  where  the  multichannel  filter  systems  without 
gain  fluctuation  were  separating  noise  and  signal  strongly  on  the  basis  of 
gain  inequalities,  the  S/N  improvement  (measured  noise  and  artificial 
infinite  velocity  spike  signal)  was  considerably  better  than  the  counterpart 
systems.  At  frequencies  where  the  systems  were  operating  on  the  basis 
of  velocity  differences  and  the  noise  generally  was  coherent  across  the 
array,  the  S/N  improvements  of  the  counterpart  systems  generally  are 
similar.  At  certain  higher  frequencies  (f  ^2  cps)  the  systems  with  gain 
fluctuation  give  notably  poorer  S/N  improvement  than  the  systems  without 
gain  fluctuation.  The  reason  for  this  is  not  apparent. 

Figure  IV-11  shows  the  K-space  response  of  IP  22  and  IP  13 
at  0.  5  cps.  IP  13  shows  very  little  noise  rejection  at  0.  5  cps,  whereas, 

IP  22  rejects  noise  strongly  in  a  velocity  band  around  3  km/ sec.  These 
responses  indicate  strongly  that  the  system  with  gain  fluctuation  <IP  22) 
was  forced  to  design  primarily  on  the  basis  of  velocity  differences  in  the 
signal  and  noise. 

Figure  IV-12  shows  the  signal -to -noise  improvement  for  the 
noise  used  in  the  MCF  system's  design  of  IP  13  and  IP  22.  (It  should  be 
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noted  that  these  S/N  improvements  were  obtained  from  spectral  estimates 
of  different  resolution.)  At  low  frequencies  the  S/N  improvement  obtained 
from  IP  13  is  much  better  (10  db)  than  that  obtained  from  IP  22.  In  view  of 
the  K -space  responses  of  the  two  systems  at  0,  5  cps  (Figure  IV-11),  the 
S/N  improvement  obtained  from  IP  13  is  false.  To  a  large  extent,  IP  13 
suppressed  the  noise  on  the  basis  of  gain  differences.  A  multichannel 
filter  system  which  would  pass  infinite  velocity-  signals  as  actually  recorded 
could  give  S/N  improvement  at  low  frequencies  on  the  order  of  that 
obtained  from  IP  22. 

Figure  IV-13  shows  the  K-space  response  of  IP  13  and  IP  22  at 
1.0  cps.  Both  systems  show  attenuation  of  plane  wave  energy  centered 
around  3  km/ sec,  but  the  system  with  gain  fluctuation  attenuates  much  more 
strongly.  This  indicates  that  while  IP  13  was  still  affected  by  gain  inequality 
problems,  it  was  obtaining  its  noise  suppression  at  1.0  cps  somewhat  more 
on  the  basis  of  spatial  organization  and  less  on  the  basis  of  gain  inequalities. 
This  also  is  indicated  by  the  frequency  response  of  the  individual  filters 
(Figures  IV-2  and  -3).  The  individual  channel  gains  are  considerably  lower 
at  1  cpf  tnan  at  0.  5  cps. 

The  S/N  improvement  at  1.  0  cps  is  only  slightly  better  for 
IP  13  than  for  IP  22.  IP  22  was  forced  to  rely  upon  velocity  differences 
in  the  signal  and  noise  and  did  a  better  job  than  IP  13  of  suppressing  noise 
on  that  basis.  IP  13,  however,  was  able  to  use  gain  differences  in  addition 
to  the  spatial  organization  of  the  noise  with  the  net  effect  beiflgjjfcfr»IP  13 
gave  a  little  better  (falsely  so)  S/N  improvement.  * 

Figure  IV-14  shows  the  K-space  power  response  of  IP  24 
and  IP  34  at  0.  3  cps.  IP  34  (system  with  gain  fluctuation)  shows  a  stronger 
reject  band  at  a  bit  lower  velocity  than  does  IP  24.  This  strongly  suggests 
that  part  of  the  noise  rejection  obtained  from  IP  24  resulted  from  using  gain 
inequalities.  This  also  tends  to  confirm  the  observation  that  high  gains  on 
individual  channels  can  mean  that  the  MCF  system  is  using  gain  differences  to 
separate  artificial  signal  and  noise. 

The  S/N  improvements  fcr  the  MCF  systems  IP  23  through 
IP  2o  and  IP  33  through  IP  36  are  shown  in  Figures  IV-1S  and  IV-16. 

Figure  IV-15  shows  the  S/N  improvement  relative  to  a  single  seismometer 
obtained  from  noise  sample  A  5  that  was  processed  by  IP  23  through  IP  26 
and  their  counterpart  systems  IP  33  through  IP  36.  Figure  IV-16  shows 
the  S/N  improvement  relative  to  a  single  seismometer  obtained  from  noise 
sample  I  5  that  was  processed  by  IP  23  through  IP  26  and  IP  33  through 
IP  36. 

The  difference  in  S/N  improvement  between  IP  34  and  IP  24  is 
8  to  9  db  at  0.  3  cps  for  noise-sample  I.  A  portion  of  this  is  due  to  the 
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Figure  IV-13.  K-Plane  Responses  of  IP  13  and  22  at  1.0  CPS 


poorer  signal  response  of  IP  34  (the  signal  responses  are  not  included 
in  this  summary).  The  major  reason,  for  this  difference,  however  is 
that  IV  24  was  using  seismometer  output  inequalities  to  suppress  the  noise. 

Figures  IV- 17  nr,  IV-18  show  the  waver  umber  responses  of 
IP  23  chrough  IP  26  and  IP  33  through  IP  36  at  1 . 0  cps,  The  counterpart 
systems  are  sem  to  have  quite  similar  K-space  responses  at  ]  ,0  tps. 
Indeed,  all  the  systems  have  quite  similar  K-space  responses  at  1 . 0  cps. 

At  this  would  indicate,  the  signal-to-noise  improvements 
are  very  nearly  equal  for  all  MCF  systems  at  1 . 0  cps  (Fi  ires  IV-I5  and 
IV-16). 


The  effects  of  multichannel  filter  systems  which  have  gain 
inequalization  problems  on  actual  signals  can  be  seer  from  ar.  event  from 
Crete  recorded  at  CPO  on  March  4,  1963,  which  is  shown  in  Figures  IV-19 
and  IV-20.  Figure  IV-19  shows  the  outputs  of  IP  23  through  IP  26  and 
Figure  IV-20  shows  the  outputs  of  IP  33  through  IP  36.  The  figures  show 
the  center  seismometer  output,  the  multichannel  filter  system  output  and 
the  difference  of  these  two.  This  event  occurred  at  a  distance  of  about 
83  degrees  from  CPO;  thus,  the  signal  should  have  an  apparent  velocity 
across  the  array  of  about  22  km/sec. 

For  this  event,  the  energy  seems  to  be  concentrated  in  a 
frequency  range  around  0.5  cps.  Such  a  signal  should  be  passed  by  all  the 
multichannel  filter  systems  with  very  little  attenuation. 

The  systems  without  gain  fluctuation  which  were  designed 
using  the  lax  ?  'r  number  of  seismometers  (IP  23,  IP  24  and  IP  25)  show 
an  abnormally  low  signal  amplitude.  They  are  low  as  compared  to  their 
counterparts  with  gain  fluctuation  and  also  as  compared  to  IP  26.  While 
this  possibly  could  be  explained  on  the  basis  of  the  proper  combination  of 
seismometer  gain  outputs,  it  strongly  suggests  that  the  design  of  systems 
TP  23,  IP  24  and  IP  25  was  affected  by  gain  inequalities.  These  effects 
appear  to  be  moderate,  but  noticeable  at  about  0.5  cps. 

E.  COMPARISON  OF  MCF  SYSTEMS  WITH  A  SUMMATION  FOR  SPECIFIC 
EVENT 


In  order  to  compare  the  multichannel  filter  systems  with  the 
summation,  a  straight  summation  of  the  19  seismometer  outputs  was  made 
for  the  Crete  event.  This  is  shown  in  Figure  IV-2T. 

The  summation  does  not  reject  the  pre-event  noise  nearly  as 
v/ell  as  the  multichannel  filter  systems  which  use  a  comparable  number  of 
seismometers.  The  summation  is  passing  a  large  ;  mount  of  low-frequency 
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Figure  IV-19.  Crete  Earthquake  Signal  as  Processed  by  IP  23,  24,  25,  and  2$ 
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Figure  IV-20.  Crete  Earthquake  Signal  as  Processed  oy  IP  33,  34,  35,  and  36 


Figure  IV-21.  Summation  of  19 -Seismometer  Outputs  for  the  Crete  Event 


{£  1/3  cps)  noise  which  the  MCF  systems  reject.  Tne  superiority  of 

multichannel  filter  systems  for  rejecting  noise  at  CPO  has  been  discussed 
previously  (  Section  V) 

F.  A  COMMENT  ON  THE  EFFECTIVENESS  OF  PARTIAL  ARRAYS  AT  CPO 

The  results  of  this  sttidy  also  lead  to  some  conclusions  as  to 
the  effectiveness  of  the  partial  arrays  at  CPO.  It  previously  was  thought 
that  the  MCF  systems  IP  23.  IP  24,  IP  25.  and  IP  26  might  have  been  ^ 

affected  significantly  by  seismometer  gain  inequalities  below  1.0  cps  (p.  35)  . 
This  now  is  indicated  to  be  the  case  only  below  about  0.  5  cps.  Some  of  the 
noise  rejection  below  0.  5  cps  shown  by  IP  23.  IP  24  and  IP  25  undoubtedly 
is  false.  Indications  are  that  below  1.5  cps,  a  properly  designed  multichannel 
filter  system  using  only  the  center  seismometer  and  the  outer  ring  would 
be  about  as  effective  as  one  designed  using  the  entire  o-rray.  However, 
the  pre -event  noise  (Figures  TV-19  ano.  IV-20)  clearly  indicates  that 
suppression  of  low-frequency  (f  <  0.  5  cps)  noise  is  somewhat  a  function 
of  the  number  of  seismometers  used. 
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SECTION  V 


A  MCF  SYSTEM  DESIGNED  FOR  THE  TFO  CROSS  ARRAY 

This  is  a  summary  of  work  which  will  be  published  as  a 
special  report  umder  the  title  "An  Evaluation  of  the  Properties  of  the 
TFO  Cross  Array  Based  on  the  Evaluation  of  a  Multichannel  Filter 
System  Designed  Using  the  Entire  Cross  Array." 

.A.  DESIGN  PROCEDURE  AND  PURPOSE 

A  21 -channel  multichannel  filter  system  was  designed  in  the 
time  domain  using  measured  noise  correlations  and  an  infinite  velocity 
signal  model.  The  channels  were  the  outputs  of  the  21  vertical  seismometers 
of  the  cross  array  {Figure  V-l).  The  signal  was  to  be  extracted  at  the 
center  seismometer  (Z-21).  The  noise  correlations  used  were  frein  the 

Tf  O  long  noise  sample.  ^  Statistical  gain  fluctuation  was  incorporated  into 
the  design  equations. 

The  multichannel  filter  oyclem  was  designed  to  determine  the 
possible  usefulness  of  multichannel  filtering  to  enhance  distant  mantle  P- 
wave  signals  at  TFO.  Designing  a  system  which  has  to  preserve  only 
infinite  velocity  signals  and  then  applying  it  to  the  noise  used  in  the  design 
should  maximize  the  MCF  system's  noise -suppressing  abilities. 

Noise  suppression  obtained  from  the  multichannel  filter 
system  discussed  in  this  report  should  represent  nearly  the  maximum 
noise  suppression  obtainable  using  thp  TFO  cross  array  without  going  to 
directional  systems.  This  MCF  system,  thus,  puts  a  rough  bound  on 
noise  suppression  possible  at  TFO.  This  multichannel  filter  system  should 
net  be  expected  to  pass  all  plane  ./ave  energy  with  velocities  *  12  km/sec. 

It  should,  however,  indicate  whether  or  not  multichannel  filtering,  using  the 
TFO  cross  array,  has  a  chance  of  offering  advantages, 
of 

B.  EVALUATION  OP  MCF  SYSTEM  AND  COMPARISON  WITH  SUMMATION 

Figure  V-2  shows  the  signal-to-nuiae  improvement  obtained 
by  the  MCF  system  and  a  summation  for  an  infinite  velocity  signal  input. 

The  signal  response  of  the  MCF  system  is  so  nearly  unity  that  the  S/N 
improvement  essentially  is  a  measure  of  noise  suppression  only.  It  can 
be  seen  that  the  MCF  system  does  a  better  job  of  suppressing  the  noise 
out  to  about  1.  2  cps;  from  1. 2  to  2e  315  cps,  the  summation  and  MCF 
system  show  about  equal  ability  to  suppress  noise.  Below  0.  9  cps,  the 
MCF  system  gives  S/N  improvement  generally  3  to  5  db  greater  than  a 
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Figure  V-l.  Arrays  at  TFO 


summation.  At  about  0.2  cps,  the  MCF  system  show?  a  signal -to -noise 
improvement  about  8  db  superior  to  the  summation. 

A  study  of  the  K-space  iesponses  of  the  MCF  system  and 
the  summation  taken  with  what  is  known  about  the  K-space  distribution  of 
the  noise  power,  gives  a  good  insight  into  the  signal -to -noise  imp  overnents. 

One  of  the  striking  features  of  the  S/N  improvements 
(Figure  V-2)  is  the  relatively  great  improvement  provided  by  the  MCF 
system  in  a  narrow  band  around  U.  2  cps.  Correlation  analysis^  indicates 
that  for  the  frequency  range  of  0  to  0.  386,  the  noise  correlations  (and, 
therefore,  filter  design)  were  influenced  strongly  by  a  Rayleigh -mode 
wavetrain  propagating  in  a  direction  about  S  60°  W  with  a  velocity  of  about 
3.  4  km/sec. 

The  K -space  response  of  the  MCF  sy  stem  at  0.  25  cps 
(Figure  V-3)  indicates  that  the  system  designed  to  "tune  out"  this  Rayleigh 
wavetrain.  A  plane  wave  with  velocity  3.4  km/sec  and  direction  S  bO°W 
would  be  attenuated  by  more  than  20  db.  Even  though  the  K-space  plot 
probably  is  at  a  slightly  higher  frequency  than  the  frequency  of  the 
Rayieigh  wavetrain,  the  K-space  response  will  vary  continuously  as  a 
function  of  frequency;  Figure  V-3  should  be  very  close  to  the  system’s 
K-space  response  at  0.2  cps. 

The  summation  K-space  response  (Figure  V-4)  indicates 
that  a  summation  would  attenuate  the  above  wavetrain  jy  about  6  db. 

At  0.25  cps,  neither  the  summation  nor  the  MCF  system  is  attenuating 
any  high  velocity  (V  ^8  km/sec)  plane  waves.  The  MCF  system  does 
suppress  some  of  the  energy  in  the  velocity  range  of  4  km/sec  5  V  5 8 
km/sec.  The  summation,  howrver,  passes  with  only  moderate  attenuation 
plane  waves  with  velocities  a.:  low  as  4  km/sec. 

o  5 

The  noise  power  at  these  low  frequencies  is  indicated  J 
to  be  made  up  essentially  of  the  Rayleigh-mode  wavetrain  pad  high- 
velocity  noise. 

Around  0.25  cps,  the  major  difference  between  the  summation 
and  the  MCF  system  is  that  the  latter  suppresses  the  Rayleigh -mode 
wavetrain  (which  was  a  part  of  the  noise  used),  much  more  strongly  than 
does  the  summation. 

Since  the  same  noise  was  used  to  design  and  evaluate  the 
MCF  system,  it  could  suppress  the  noise  power  of  the  Rayleigh  wave 
exceptionally  well.  This  gave  the  substantially  greater  S/N  improvement 
at  0.2  cps.  This  MCF  system  would  not  be  expected  to  show  similar  S/N 
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improvement  at  0.  2  cps  on  another  noise  sample  which  did  not  include 
such  a  directional  Rayleigh  wa /strain. 

At  C.  5  cps  the  MCF  system  is  attenuating  a  larger  portion 
of  the  high-velocity  energy  (Figure  y-5).  Since  the  noise  at  0.  5  cps  is 
indicated  to  be  mostly  high  velocity  one  would  expect  the  S/N 

improvement  of  the  MCF  system  to  be  noticeably  better  than  the  summation 
at  0,  5  cps.  Figure  V-2  snow3  that  the  MCF  system  gives  a  S/N  improve¬ 
ment  about  4  db  superior  to  a  summation  at  0.  5  cps. 

The  K -space  plots  oi  the  MCF  system  and  the  summation 
at  1.0  cps  (Figures  V-6  and  V-4)  show  that  the  MCF  system  is  suppressing 
more  high-velocity  energy  at  this  frequency  also.  Both  the  MCF  system  and 
the  summation  are  suppressing  more  of  the  high-velocity  energy  than  they 
v/ere  at  0.  5  cps.  The  summation  is  attenuating  substantially  some  of  the 
high-velocity  plane  wave  energy  at  1.0  cps;  little  attenuation  was  obtained 
at  0. 5  cps. 


3  5 

Previous  studies  of  this  TFO  noise  ’  indicate  that  the 
noise  a '  1.0  cps  is  predominately  high  velocity  (V^8  km/sec).  Thus,  based 
on  the  K -space  responses,  one  would  expect  that  the  MCF  system  would 
give  somewhat  better  S/N  improvement  than  the  summation  at  1.0  cps. 

Also,  both  systems  should  show  better  S/N  improvement,  at  1.0  c{.s  than  at 
0. 5  cps. 


Figure  V-2  shows  that  the  actual  difference  at  1.0  cps  is  1  to 
2  db  and  that  both  systems  do  give  better  S/N  improvement  at  1.  0  cps  than 
at  0.  5  cps. 


Above  i.  2  cps  the  S/N  improvement  of  the  MCF  system  and 
the  summation  are  very  similar.  Noise  analysis  5  indicates  that  this  TFO 
noise  sample  generally  is  uncorrelated  above  1.  2  cps.  What  correlation 
there  is  above  1.2  cps  is  not  space -stationary.  Thus,  above  1.2  cps,  any 
attempt  to  explain  the  signal-to-noise  improvement  on  the  basis  of  the 
K-space  response  is  not  possible. 

If  the  noise  were  random,  the  summation  should  give  a 
signal-to-noise  improvement  of  13.2  db.  However,  for  the  noise  sample 
used  in  evaluation,  the  signal-to-noise  improvement  of  the  summation  and 
♦he  MCF  system  (Figure  V-2)  in  the  frequency  range  of  1. 5  to  2.  3  cps 
generally  is  somewhat  greater  than  13.2  db.  At  certain  frequencies,  it  is 
much  better.  In  a  narrow  band  around  1.6  cps,  the  S/N  improvement  of 
both  the  MCF  system  and  the  summation  is  about  21  dc. 

\Vhether  or  not  one  generally  could  expect  S/N  improvements 
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Figure  V-6. 


K-Plane  Response  of  Multichannel  Processor  at  1.  00  CPS 
(Contoured  in  db) 
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at  these  frequencies  in  excess  of  the  theoretical  for  random  noise  is  not 
known. 


The  MCF  system  is  better  for  suppressing  noise  below 
1. 2  cps  and  particularly  superior  in  a  band  around  0.  2  cps.  Above  1.  2 
cps  the  MCF  system  and  thn  summation  suppress  the  TFO  ambient  noise 
about  equally. 


The  exceptionally  good  results  obtained  by  the  MCF  system 
in  a  narrow  band  around  0.  2  cps  is  the  result  of  suppressing  some  strongly 
directional  Rayleigh -mode  energy  in  this  pa^ticu  ar  noise  sample.  The 
MCF  system  would  not  show  this  good  S/N  improvement  at  G,  2  cps  when 
applied  to  noise  which  contained  no  such  directional  low-velocity  energy 

The  3  to  5-db  increased  S/N  improvement  obtained  from  the 
MCF  system  below  1.0  cps  generally  is  the  result  of  better  suppression  of 
high-velocity  (V  5  8  km/sec)  noise.  The  MCF  system  was  designed  using  an 
infinite  velocity  signal  model  and,  therefore,  tried  to  eliminate  arv 
coherent  energy  of  lower  velocity.  The  MCF  system  was  able  to  cut  more 
sharply  in  K-space  near  V  =  “  than  was  the  summation.  This  effect  is 
significant  at  low  frequenciec  and  resulted  in  greater  S/N  improvement  for 
the  MCF  system  as  compared  with  the  summation. 

In  the  course  of  this  study,  several  tests  of  the  stability  of 
the  summation  and  the  multichannel  filter  system  to  perturbation  of 
seismometer  gains  were  made.  These  tests  indicate  that  both  are  quite 
stable  when  applied  to  noise  for  which  the  seismometer  outputs  have  been 
perturbed  moderately  (25  percent).  The  details  of  these  results  are  omitted 
here,  but  will  be  published  in  the  forthcoming  special  report  covering  this 
study. 

C.  A  DISCUSSION  OF  ISOTROPIC  PROCESSING  AT  TFO 


Because  of  the  K-space  distribution  of  the  TFO  noise 
(largely  high  velocity  with  some  surface-wave  energy),  isotropic  pro¬ 
cessing  is  not  very  effective.  Any  processor  which  must  maintain  a  good 
signal  pass  response  for  energy  propagating  with  a  velocity  of  1  2  to  24 
km/sec  cannot  suppress  the  high-velocity  noise.  An  isotropic  processor 
could  be  expected  to  suppress  only  the  surface  -wave  energy  indicated  tc 
be  present  in  at  least  part  of  the  TFO  noise  and  the  random  noise. 
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Except  in  a  narrow  band  around  0,  2  cps,  the  better  low- 
frequency  (f  <1.  0  cps)  S/N  improvement  obtained  from  the  MCF  system,  as 
opposed  to  the  summation,  resulted  from  rejection  of  high-velocity  noise. 
This  was  possible  because  the  multichannel  filter  system  was  "directional" 
(i.e. ,  an  infinite  velocity  signal  model).  Such  noise  rejection  would  not 
be  possible  with  any  system  which  maintained  a  suitable  signal  response 
in  the  range  of  V  s  12  km/sec. 

Isotropic  multichannel  filter  processing  at  TFO  would  be 
effected  little  by  the  array  configuration.  Any  reasonable  configuration 
could  be  used  to  suppress  the  surface-wave  energy.  The  suppression  of 
the  random  noise  is  a  function  of  the  number  of  seismometers,  not  their 
location. 


On  the  basis  of  this  study  and  the  various  TFO  noise  studies 
referred  to,  it  is  implied  that  the  effectiveness  of  isotropic  multichannel 
filtering  at  TFO  should  not  be  expected  to  be  affected  much  by  array 
geometry.  Designing  isotropic  multichannel  filter  systems  for  partial 
arrays  at  TFO  would  give  little  useful  information. 

D.  A  DISCUSSION  OF  DIRECTIONAL  MULTICHANNEL  FILTERING  AT  TFO 

The  multichannel  filter  system  discussed  in  this  report  is  a 
"directional"  system  (i.e.  ,  the  system  was  designed  to  pass  an  infinite 
velocity  signal  only).  This  -'directional"  system  was  superior  to  the 
summation  at  frequencies  below  about  1.2  cps.  At  these  lower  frequencies 
the  multichannel  filter  system  was  superior  in  that  it  suppressed  a  portion 
of  the  high-velocity  noise  better  than  the  summation  and  it  did  a  better  job  of 
suppressing  some  surface-mode  energy  present  in  the  noise  sample  pro¬ 
cessed.  Below  about  1.  0  cps,  the  MCF  system  delivers  about  3-to  5-db 
better  suppression  of  the  mantle  P-wave  noise  than  the  summation. 

4  5  6 

TFO  noise  analysis  ’  ’  indicates  that  the  low-frequency 
ambient  noise  is  mostly  mantle  P-wave  noise.  If  this  noi3e  can  be  roughly 
modeled  as  a  disc  of  uniform  power  in  K -space  with  a  minimum  velocity  of 
about  8  km/ sec,  one  would  expect  any  directional  MCF  system  to  suppress 
the  mantle  P-wave  noise  at  least  as  well  as  a  system  designed  using  an 
infinite  velocity  signal  model.  If  a  multichannel  filter  system  were  designed 
to  pass  a  directional  signal  not  in  the  center  (in  K-spacO  of  this  high-velocity 
noise  field,  it  might  be  considerably  more  effective  than  the  system  designed 
using  an  infinite  velocity  signal  model.  If  the  directional  signal  model  were 
nearer  the  edge  of  the  noise  power  d  3C  in  K -space,  the  signal  pass  region  of 
the  processor  could  lie  partially  outside  the  noise  field  at  frequencies  of 
interest.  Also,  such  a  directional  multichannel  filter  should  develop  an 
asymmetric  K -plane  response  which  would  further  attenuate  the  high-velocity 
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Similarly,  a  beam -steering  would  be  expected  to  suppress 
the  mantle  P-wave  noise  as  well  as  or  better  than  the  straight  summation, 
but  would,  of  course,  maintain  a  symmetric  K-space  response  around  the 
point  of  unit  response. 

Any  low-velocity  energy  could  be  suppressed  effectively  by 
a  properly  designed  directional  multichannel  filter  system.  The  effectiveness 
of  the  time-shifted  sum  for  suppressing  surface-mode  energy  would  depend 
upon  the  K-spacc  response  of  the  particular  Learn -steering,  but  would  be 
expected  to  be  inferior  to  a  MCF  system. 

Thus,  for  directional  processing  using  the  cross  array  at 
TFO,  one  would  expect  that 

a  Below  about  1.0  cps,  a  properly  designed  multichannel 
filter  system  shoul  i  be  at  least  3  to  5  db  better  (in  S/N 
improvement)  than  a  beam -steering. 

•  When  applied  to  noise  containing  coherent  low -velocity 
energy  (surface -mode  energy),  the  nun iichannel  filter 
system  might  give  much  greater  S/N  improvement  than 
a  beam -steering.  A  study  of  directional  multichannel 
filtering  using  the  TFO  cross  array  is  now  in  progress. 
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SECTION  VI 


A  SUMMARY  OF  SIGNAL  DETECTION  PROCESSING  RESULTS 

The  assumptions,  derivations  of  equations,  and  an 
illustration  of  the  method  of  signal  detection  by  probabilistic  processing 
were  presented  in  Special  Report  No.  8,  The  illustration  used  artificially 
generated  data  to  satisfy  the  assumptions.  The  expected  results  were 
determined  theoretically  and  compared  favorably  v/ith  the  actual,  processor 
output. 


The  basic  assumption  is  that  noise  and  signal  data  are 
independent  multivariate  normal  with  zero  means  and  known  covariance 
matrices  given  by  fi  and  0 respectively.  The  elements  of  both  covariance 

matrices  were  generated  by  an  existing  program  from  the  following  models 
of  signal  and  noise: 

e  The  3 -dimensional  noise  spectrum  in  (f,  k  ,  k  )  space 

x  y 

is  assummed  to  be  zero,  except  for  points  corresponding 

to  velocities  between  3  km/sec  and  4  km/sec.  The 

spectrum  is  uniform  for  values  of  f,  k  and  k  in  the 

x  y 

velocity  interval. 

•  The  3 -dimensional  signal  spectrum  is  assumed 

to  be  zero,  except  for  values  of  f,  k  and  corresponding 

to  velocities  greater  than  8  km/sec  and  constant  in  the 
high-velocity  region.  The  zero-lag  autocorrelation  values 
were  set  equal  and  the  spectra  remained  "white  in  frequency 
for  both  signal  and  noise  in  the  previous  report. 

This  section  will  contain  some  graphical  information  related 
to  Special  Report  No.  8  and  some  results  showing  the  theoretical  effect  of 
varying  the  signal -to-noise  ratio.  The  results  of  applying  the  quadratic 
processor  to  three  sections  of  10 -channel  C00  data  containing  teleseismic 
events  will  then  be  discussed. 

A.  GRAPHS  RELATED  TO  SPECIAL  REPORT  NO  8 

The  simple  models  {4 -channel,  annular -ring,  3  to  4  km/sec 
noise  and  8-km/sec  solid  disk  signal)  were  used  to  generate  correlation 
data  and  to  produce  from  this  data  the  matrix  constants  of  the  quadratic 
form.  The  distributions  of  the  quadratic  processor,  under  the  two 
hypotheses  of  I  =  noise  and  II  =  signal  +  noise,  were  shown  as 
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(6.1) 


P  (x)  =  6  x  ** 

r  '  i  vj 


PH(I!!  4  6II  x  V 


6  =  o7zn 

2  2 

v  =  24/0 


(6.2) 


and  the  values  4,  o'  are  the  mean  and  variance  of  the  processor  output 

P(x),  and  x^  is  a  chi -squared  random  variable  with  V  degrees  of  freedom. 

V 

Plots  of  the  cummulative  distributions  of  Equation  (6.  1)  are  given  in 
Figure  VI - 1  for  10  and  25  lags.  The  vertical  lines  at  C{10)  =  25.  3  and  C(25) 
=  6?.  5  are  the  equi -probability  critical  values. 

B.  EFFECT  OF  VARYING  SIGNAL -TO -NOISE  RATIO 


There  are  actually  two  questions  involving  signal -to -noise 


ratios; 


1.  Hew  will  P(x)  designed  for  S/N  =  1  behave  for  data 


x  «  n2  =  N  +  aS? 

2.  How  does  the  separation  of  I  and  11(a)  depend  on  a; 
here  P(x)  is  dependent  on  a. 

Considering  the  first  problem:  since 
P(x)  =  yT(.\-  I)y 

where  x  -  My  and  x  €  O^T  (a)  =  N  +  aS,  then  y  has  a  multivariate  normal 
distribution  with  mean  zero  and  covariance 

M"1  T  =  M_1  (a  [s  +  n|  +  (1  -  a)N)  (M_1)T 


=  al  +  (1  -  a)  A  =  F 
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Figure  VI -1.  Cummulative  Distribution  Plot  of  4-C-hannel  10  and  25  Lags 
for  Noise  (I)  and  Signal  Plus  Noise  (II) 
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Now,  transform  y  to  z  by  y  =  F  zso  that 

f,  =E,z’r--F-1/2FF-1/2=I 

7. 


and 


P  (y)  =  zT  (F1/2A  F1/2  -  F)z 


Tl'us,  P[y)  has  the  same  distribution  as 


P  2 

I  n.  x,  (i) 

i  =  l  1  ' 

with 

a  |a  +  (1  ~  ?)  X\y  \  "  l) 

where  X .  is  the  diagonal  w-leir.ent?  o£  A  and  Xj  =  1/X..  Thus,  the  mean  and 
variance  when  x  z  =  N  +  aS  are 

(J*  =  aMjj  +  (1  -  a)  iij 


and 

°z *  =  a2  a2  +  (1  -  a)  a2  +  4a(  1  -  a)  -  Mj) 

The  behavior  of  the  quadratic  processor  when  applied  to 
data  of  varying  S/N  ratio  can  be  described  completely  in  terms  of  the 
beha’dor  on  noise  alone  anl  signal -plus -noise  with  the  ratio  1, 

The  second  question  does  not  have  as  complete  an  answer. 
Suppose  we  adjust  the  S/N  ratio,  thereby  obtaining  iXj  =  N  +•  aS,  and  consider 

this  8.  new  design  problem.  It  is  desirable  to  determine  the  behavior  of  the 
processor  for  this  new  problem  in  terms  of  the  results  of  the  old  problem. 
Just  as  before,  if  X'":  satisfies 


Then,  the  behavior  of  the  quadratic  processor  is  specified  completely 
by  the  quantities 
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and 


E  X  E  X  “  when  x  s  II  (a) 


❖  .  1  * .  2 

E  X,  E  X  when  x  e  I 

i  i 


-  X*  ^1  =  0  is  the  same  equation  in  X*  as  IfL,  -  X  0^ 


Now  In'1 
which  can  be  written 


=  0 


* 

N  +  aS-X  N  =0 


SN  -  {X"  -  .1)  /  al  |  = 


so  that 


and  thus 


a.  -1  =  {  X  -1) /a 

i  i 


u‘C  =  ft  V 


II 


2*  2  n2 

c  =  a  o 


II 


so  that  the  behavior  n  P(x)  as  function  of  a  can  be  described  when  x  e  11(a). 
Thf>  behavior  of  P(x),  whenxe  I,  has  not  been  solved  as  a  simple  function 
of  a. 

The  curves  drawn  in  Figure  VI -2  illustrate  the  above 
equations  for  the  lower  (Zb  lags)  example  of  Figure  VI -1.  The  curves  in 
the  upper  portion  of  Figure  VI -2  show  the  behavior  of  P(x)  when  x  «  aS  +  N 
for  a  =  1/4,  1  /2,  3/4,  and  the  processor  P(x)  is  designed  for  a  =  1.  The 
lower  portion  of  Figure  VI -2  shows  the  behavior  of  P(x)  for  a  =  1/4,  1/2, 

3/4  when  the  processor  is  designed  for  these  values. 

C.  PROCESSING  OF  ACTUAL  DATA 

Three  sections  of  10-channel  CPO  data  were  selected 
containing  weak  teleseismic  signals  with  the  following  descriptions: 
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Figure  VI -3  shows  the  location  of  the  seismometers  for  the  10  channels  of 
data  used  as  input  to  the  quadratic  processor.  The  noise  model  was  assumed 
to  be  3  to  4  km/sec  annular  ring  and  an  8-km/sec  solid  disk  signal  model 
was  used.  Ten  percent  white  noise  was  added  to  the  noise  correlation 
matrix. 

The  output  of  the  quadratic  processor  for  the  three  sections  of 
data  is  presented  in  Figures  VI-4,  VI -5  and  VI  -6.  The  traces  shown  are  the 
quadratic  process,  Wiener  square  and  sum  (over  8  points),  and  the  actual 
seismometer  traces.  The  number  of  lags  used  for  both  the  quadratic  pro¬ 
cessor  and  the  Wiener  output  was  17  with  no  smoothing  for  the  quadratic 
processor. 

Tne  theoretical  critical  level  is  drawn  in  for  the  quadratic 
detec  ton  processor  on  all  three  sections  (Figures  \T-4,  VI -5  and  VI-6). 

An  error  apparently  was  introduced  at  some  stage  of  the  calculation  of  the 
critical  level  for  the  second  section.  The  other  two  critical  levels  are 
believed  to  be  correct  and  typical  of  what  might  be  expected  in  additional 
noise  samples.  Two  of  the  teleseismic  everts  did  not  produce  a  detector 
output  significantly  greater  than  the  noise  output.  The  third  section  of 
data  shows  a  significant  detector  output  for  the  teleseism.  The  signals 
’ater  in  this  section  (Figure  VI -6)  are  the  arriving  phases  of  a  quarry  blast 
which  also  produced  a  significant  value  of  the  detector  output. 

The  results  for  the  detection  processor  show  sufficient  promise 
to  continue  this  study  using  the  full  19-channel  CPO  array  with  ring  data 
and  actual  -measured  CPO  noise  correlations. 

The  theoretical  cummulative  probability  curves  for  the  models 
of  signal  and  noise  used  are  shown  in  Figure  VI  -7.  Some  cummulative 
points  for  the  actual  data  were  determined  from  the  first  500  points  of  the 
first  section  and  were  plotted  in  Figure  VI -7  for  comparison. 
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Figure  VI -4C  Quadratic  Processor,  Wiener  Trac;  and  Raw  Data  for  Section 
1  of  10 -Channel  CPO  D  ta 
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Figure  VI '7,  Theoretical  and  Actual  Cummulative  Probability  Curves 
for  the  10 -Channel,  17 -Lag  CPO  Data 
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SEC 'HON  VII 


LARGE  APERTURE  SEISMIC  ARRAY  DETECTION  FILTERS 
A.  SUMMARY 


Several  sets  of  multichannel  filters  were  developed  in 
order  t  >  investigate  the  nuclear  detection  capabilities  of  the  21  -element* 
Large  Aperture  Seismic  Array  (LASA)  located  in  Montana.  Each  filter 
was  developed  to  pass  a  signal  from  a  known  cr  possible  nuclear  te^t 
site  and  to  reject  noise  events  from  a  known  arc;a  of  earthquake  activity. 

The  important  factor  of  this  technique  of  processing  the  LA3A  is  that  the 
signal  area  is  allowed  to  be  monitored  with  a  single  multichannel  filter, 
while  many  beam  steers  would  be  required  to  monitor  the  same  area. 

In  developing  the  LASA  multichannel  filters,  two  noise 
models  and  three  signal  models  were  chosen.  The  first  of  the  noise 
models  included  the  portion  of  world  in  £-space  within  a  velocity  greater 
than  or  equal  to  12  km/sec,  while  the  second  noise  model  encompassed 
the  major  earthquake  regions  of  the  earth.  The  sigtal  models  chosen 
were  Russia,  China  and  Novaya  Zemlya.  These  models  encompassed  some 
of  the  possible  or  known  nuclear  test  sites.  Four  of  the  multichannel  filters 
designed  which  used  various  combinations  of  the  signil  and  noise  models 
are  presented  in  this  section. 

The  results  of  this  study  indicated  that  t  is  possible  to 
monitor,  within  a  0-  to  6-db  passband,  a  reasonable  signal  area  with  a 
single  filter  while  retaining  the  ability  to  reject  3  to  3G  db  or  more 
known  earthquake  areas  if  some  tr  de-off  in  random  noise  response  is 
allowed.  This  implies  that  any  suspected  nuclear  test  crea  may  be 
mouicored  with  one  multichannel  filter  instead  of  many  beam  steers. 

Such  a  system  designed  on  this  basis  not  only  would  simplify  considerably 
the  detection  problem  as  related  to  LASA,  but  also  would  result  in  a 
significant  monetary  savings  when  compared  with  the  beam -steer  technique. 

Based  on  the  results  of  this  approach  to  monitoring  specific 
geographic  areas,  additional  research  for  improved  filter  performance  is 
warranted.  Specific  items  which  should  be  investigated  a3’e  the  following: 


e>  Use  of  additional  senerrs  tc  5"iprove  the  wide 

CT-band  filter  performance.  This  can  be  accomplished 
in  one  of  two  ways;  first,  by  considering  the  subarrays 
as  more  than  one  eleo  it  (particularly  those  located 
close  to  the  center  of  the  array)  or,  second,  by  adding 
additional  subarrays  to  the  LASA.  This  latter  approach 
has  been  suggested  previously  on  the  basis  of  an  array 
configuration  study.* 

t  Determination  of  the  optimum  signal -to -noise  ratio  and 
the  random  (in  Pf-space)  noise  level  to  be  used  in  filter 
development. 

•  Development  of  more  refined  theoretical  coherent 
noise  and  signal  models  which  wou’d  insure  more 
accurate  lc-space  area  coverage. 

•  Research  to  establish  a  more  definite  relation  between 
signal  area  coverage  and  noise  rejection  capability. 

Such  information  would  prove  beneficial  in  setting  an 
upper  limit  on  effective  (0-  to  3-db)  signal  area  coverage 
in  k-space  while  still  allowing  the  filter  to  possess 
specified  coherent  noise  response  properties. 

B.  FORMULATION  OF  NOISE  AND  SIGNAL  MODELS 

In  order  to  model  accurately  the  defined  noise  and  signal 
areas  used  in  this  study,  the  geographic  areas  of  the  world,  as  centered 
on  the  LASA,  were  transformed  to  f-k  space  to  provide  the  necessary 
ic-space  pattern.  Figure  VII -1  is  this  I?-space  map.  In  this  figure,  the 
core  shadow  region  is  the  area  from  which  no  direct  P-waves  can  propagate 
and  still  arrive  at  the  LASA.  The  small  circles  outline  the  Cir  cum -Pacific 
belt  which  is  the  source  area  for  the  majority  ox  all  recorded  earthquakes 
and  consequently  is  a  prime  rejection  area  for  filters  to  be  used  in  the 
detection  problem.  The  small  triangles  represent  known  underground  test 
sites. 


To  develop  a  model  to  fit  the  noise  and  signal  areas  of 
interest,  several  different  sized  isotropic  disks  representing  coherent, 
noise  power  spectra  were  developed  using  Bessel  functions  (Table  VII -1). 
The  equation  for  the  coherent  power  spectrum  M..(f)  of  an  annulus  area  is: 


*  Texa3  Instruments,  1965,  Large  aperture  seismic  array  final 
specification  report:  AFT  AC,  Aug.,  p,  VI-1. 
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Figure  VII- 1.  K-3pace  Mapping  of  the  World  for  P-Wave  at  1.0  CPS 
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where 

M(£)  ■-  desired  spectral  weighting 
=  first-order  Bessel  function 
x^  =  sensor  coordinate  of  th  i  sensor  (subarray) 

V  =  minimum  velocity 

=  maximum  velocity 

The  coherent  power  spectrum  of  a  disk  region  is  formed  by 

setting  V  =  or> , 


In  the  cases  where  the  signal  or  noise  areas  of  interest 
were  not  isotropic  (i.  e. ,  centered  about  £=0),  the  area  was  approximated 
by  shifting  in  £ -space  the  appropriate  models  in  Table  VII-1.  This 
shift  can  be  compared  to  a  beam  steejr  or  time  shift  of  the  coherent  noise 
(or  signal)  power  M.  .(f).  The  resulting  shifted  coherent  noise  or  signal 

i  ^ 

power  M..(f)  is  given  by: 


Ml.(f)  ^  M..(f)ei2T1k*  Xij 
ij  ij 


where 


x. . 
ij 


and 

k  =  vector  wavenumber  location  of  the  center  of  the  shifted 
model. 


’’“Texas  Instrument,  .  196 1:  Final  Rpt. ,  AFTAC  Project  VT/077,  p.  108. 


VII-4 


Table  VII- 1 


SIZE  OF  DISKS  FCk  ISOTROPIC  NOISE  AND  SIGNAL  MODELS 
(all  Rj  =  0  unless  otherwise  specified) 

Model  R  in  Velocity  Values 


A 

B 

C 

D 


12.  00  km/sec 

163,  50  km/ sec 

400.  00  km/ sec 

21.  92  km/ sec 
25.  ?2  km/ sec  R 

4 


Using  this  technique,  complex  signal  and  noise  models  may  l»e  approxi¬ 
mated  by  shifting  many  disks  of  appropriate  size  to  cover  the  desired 
area  and  then  stacking  these  crosspower  spectra,  term  by  term,  to 
form  a  single  complex  crosspower  spectrum  matrix,  approximately 
representative  of  the  desired  model. 

1,  Noise  Models 

The  first  noise  model  used,  Model  I,  was  chosen  as  simple 
as  possible  '•o  test  the  merits  of  the  system.  This  model  *.&  composed 
of  one  Model  A  disk  (Table  VII-1),  corresponding  to  a  minimum  velocity 
of  12  km/sec,  centered  at  the  origin  of  the  £-space  map  and  is  shown  in 
Figure  VII -2. 

The  second  and  more  complicated  noise  model,  Model  II, 
also  is  shown  in  Figure  VII-2.  This  model  consists  of  14  small  disks 
(Model  C),  which  were  shifted  from  the  origin  to  cover  the  Aleutian  Island 
chain,  and  one  annulus  region  (Model  D)  centered  at  the  origin.  This 
composite  model  encompasses  rr^st  of  the  Cir  cum -Pacific  belt  and  the 
southern  coast  of  Alaska. 

It  should  be  noted  that  any  earthquake  from  the  southern 
half  of  the  map  will  not  be  confused  with  a  nuclear  event  (assumed  to 
originate  from,  the  northern  half  of  the  map),  since  these  two  events  may 
be  distinguished  by  the  direction  of  propagation  of  the  waves. 
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2,  Signal  Models 


Signal  models  were  developed  to  cover  three  of  the  known 
and  possible  nuclear  test  sites.  These  models  were  developed  using  the 
same  method  and  disks  which  were  used  to  generate  the  noise  models. 


The  Jirat  signal  model  developed  consisted  of  one  Model  C, 
that  was  shifted  to  encompass  the  Russian  test  site  on  Novaya  2  ^rnlya. 

This  model  which  is  shown  in  Figure  VII -3  gave  the  best  results  by  far, 
as  will  be  shown  in  the  following  parts  of  this  section. 

The  second  signal  model  was  developed  to  encompass  all  of 
Ruscia.  This  signal  (Figure  VII-3)  simulates  Russia  with  seven  Model  W 
disks  and  nine  Model  C  disks  (Table  VII-1). 

The  third  and  final  signal  model  (Figure  VIJ.-3)  uses  seven 
Model  C  disks  (Table  VI I  - 1 )  to  simulate  China. 

c,  PRESENTATION  OF  MULTICHANNEL  FILTERS 

A  total  of  four  multichannel  filters  were  developed  using 
\arious  combinations  of  the  noise  and  signal  models  outlined  in  the 
previous  subsection.  The  details  of  these  four  filters  are  outlined  in  Table 
VII-2.  This  subsection  presents  the  wavenumber  response  at  frequencies 
of  0.  5,  1.  0,  2.  0,  and  3.  0  cps  of  these  four  filters  and  also  the  response 
of  the  filters  to  random  noise.  The  straight-sum  response  of  the  LASA 
has  been  included  (Figure  VII-4)  for  comparison  with  the  filter  results. 

Table  VII -2 
LASA  FILTERS 


Filter 

Designation 

Noise 

Model 

Signal  Model 

Signal  -to -Noise 
Ratios 

As sumed  Random 
Noise  Level 

LASA  MCF-1 

I 

Novaya  Zemlya 

8.0 

1% 

LASA  MCF -2 

II 

Novaya  Zemlya 

8.0 

5% 

LASA  MCF -3 

II 

Russia 

3.0 

5% 

LASA  MCF -4 

II 

China 

8.0 

5% 

Figure  VII-3.  Signal  .Models 


Figure  VII-4.  LASA  Straight-Sum  Response 


1,  Wavenumber  Response 


LASA  Mf*F -1  was  developed  using  he  simplest  noise 
and  signal  models,  Noise  Model  I  and  Novaya  Zernlya  Signal  Model,  to 
test  the  merits  of  the  technique,  since  these  particular  models  were  the 
least  complex  and  easiest  models  to  construct  which  would  adequately 
encompass  a  specified  signal  area  and  the  noise  areas.  The  response 
of  this  filter  is  shewn  in  Figure  VII -5  at  0.  5,  1  0,  2.0,  and  3,0  cps. 
Because  of  the  model  simplicity,  the  filter  It  able  to  strongly  reject 
the  noise  area.  This  is  best  seen  at  0.  5  cps  where  the  entire  noise 
model  is  either  12  or  IS  db  down  from  the  signal.  As  frequency  increasef, 
the  noise  rejection  deteriorates  as  evidenced  by  the  large  number  of  -6  db 
levels  at  3.0  cps.  This  property  will  be  discussed  in  detail  in  subsection 
D. 


The  second  set  of  responses  shown  in  Figure  VII-6  are 
for  LASA  MCF-2#  which  was  designed  using  Noise  Model  II  and  the 
Novaya  Zernlya  Signal  Model.  The  results  obtained  from  mis  set  of 
responses  again  are  excellent  for  the  lower  frequencies  of  0.  5  and  1.0  cp  3 
and  indicate  the  same  gradual  decrease  in  effective  results  as  the  frequency 
increases. 


The  responses  of  LASA  MCF-3  which  was  designed  using 
Noise  Model  II  and  the  Russia  Signal  Model  are  shown  in  Figure  VlI-7. 
This  filter  demonstrates  results  similar  to  the  preceeding  two  filter 
responses.  The  filter  performs  well  at  the  lower  frequencies,  but 
degenerates  as  frequency  increases.  Of  particular  note  are  the  excellent 
results  at  0.5  and  1.0  cps  where  Russia  is  at  0  and  -3  db  and  the  noise 
model  is  at  either  -12  or  -18  db  or  greater. 

The  responses  of  LASA  MCF-4  which  was  designed  on 
Noise  Model  II  and  the  China  Signal  Model  are  shown  in  Figure  VII -8.  As 
in  previous  cases,  the  results  are  very  good  for  the  lower  frequencies, 
with  the  same  decrease  in  results  as  the  frequency  increases.  It  should 
be  noted  that  the  results  for  China  are  better  than  those  for  Russia,  but 
are  not  as  good  as  those  for  Novaya  Zernlya. 


*The  shaded  areas  show  the  areas  of  prime  interest  in  this  and  the  remaining 
sets  cf  filter  responses. 
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2.  Random  Noise  Respenae 


Tne  random  noise  response  for  LASA  MCF-1  through 
MCF-4  has  been  computed  and  i*.  shown  in  Figure  VII-9.  Briefly,  the 
random  noise  response  of  a  multichannel  filter  is  given  by: 

N  (f)  N 

Niff  =  ^  I 


where 


N  =  number  of  sensors 

H.(£)  =  complex  transfer  function  of  the  filter 

N^(f)  =  input  noise  power  of  the  reference  sensor 

N^(f)  =  output  noise  power 

Knowledge  of  the  random  noise  response  of  a  multichannel 
filter  is  an  important  factor  for  consideration  when  evaluating  a  particular 
filter,  since  this  information  is  indicative  of  the  filter  performance  in  ar. 
actual  noise  field  containing  random  noise.  Comparison  of  the  random  noise 
response  should  be  made  in  terms  of  the  expected  response  for  straight- 
summation  processing  which  yields  a  noise  power  improvement  of  N  if 
the  noise  field  is  completely  random  and  is  the  upper  limit  of  improvement 
obtainable  if  the  processing  system  has  0-db  response  for  the  specified 
class  of  signal  over  the  wideband  frequency  response. 

In  line  with  this  discussion,  it  should  be  noted  that  LASA 
MCF-1  indicates  noise  rejection  greater  than  N  (-13.22  db  for  21  sensors) 
at  frequencies  greater  than  1.0  cps.  This  is  due  to  the  filter  response  to 
signals  from  Novaya  Zemlya  (the  specified  signal)  which,  above  1.0  cps, 
indicates  rejection  on  the  order  of  3  to  6  db  (Figure  VII-5).  In  other 
words,  MCF-1  is  actually  frequency  filtering  as  well  as  velocity  filtering 
which  causes  greater  than  N  rejection.  This  filter  property  probably 
resulted  because  of  the  overlapping  in  k-space  of  the  signal  and  noise 
models  and  probably  can  be  avoided  in  future  filters  by  specifying  a  larger 
signal-to-noise  ratio. 

In  order  to  lend  physical  interpretation  to  the  information 
presented  in  Figure  VII-9,  an  estimate  was  made  of  the  subarray  random 
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Figure  VII-9.  MCF  Response  to  Random  Noise 
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no’se  level  ana  the  quantity  computed.  This  information  along  with 

rT 

i 

the  estimated  random  ioise  level  is  shown  in  Figure  VII -10  and 
represents  the  anticipated  random  noise  output  for  MCF-1  through  -4 
and  for  straight  summation  relative  to  1, 0  m  of  ground  motion/cps 
at  1. 0  cps. 


The  estimated  input  random  noise  level  is  based  upon  the 
straight  summation  output  for  10  elements  of  the  Angela  subarray  which 
was  assumed  average  for  each  of  the  21  subarrays.  ThJ  j  level  was 
formed  from  power -density  spectra  and  prediction  filtering  results 
presented  in  the  report  "Large  Aperture  Seismic  Array  Final 
Specifications  Report.. "  * 

D.  DISCUSSION  OF  RESULTS 

The  wavenumber  responses  for  LASA  MCF-1  through 
-4  indicate  that  each  filter  is  passing,  with  0-to  6-db  attenuation,  the 
desired  signal  area  and  is  attenuating  the  specified  noise  regions  3  to 
30  db  and  greater.  ~n  analysing  each  of  the  wavenumber  responses, 
the  fact  tha1:  there  is  a  definite  trade-off  between  signal  Ic-space  area 
coverage  and  filter  noise  rejection  capability  is  evident.  Analysis  of  the 
random  noise  responses  probably  gives  the  best  demonstration  of  this 
trade-off.  (I.  should  be  remembered  that  this  approach  is  limited  since 
it  does  not  show  the  rejection  capability  for  the  specified  noise  regions 
but  only  general  filter  trends).  The  Novaya  Zemlya  filters  (LASA  MCF-1 
and  -2),  which  have  the  smallest  defined  signal  area  of  the  filters  being 
presented  indicate  in  Figure  VII-9  the  most  desirable  random  noise 
responses,  when  compared  with  the  other  filters.  As  the  signal  I?-space 
area  coverage  increases,  the  random  noise  response  degenerates, 
as  evidenced  by  the  China  filter  (LASA  MCF-4),  which  shows  results 
midway  between  the  Russia  filter  (LA3A  MCF-3)  and  the  Novaya  Zemlya 
filters. 


The  trade-off  between  signal  area  and  filter  noise  rejection 
can  also  be  seen  within  a  single  filter  as  frequency  increases,  since  the 
signal  (or  noise)  area  which  is  normally  defined  in  terms  of  constant 
velocity,  is  proportional  to  f2.  As  an  example,  LASA  MCF-3  (Russia 
Signal  Model,  Figure  VII -7)  shows  that  at  0.5  cps  the  noise  areas  generally 
are  attenuated  12  to  18  db  while  at  3.0  cps  this  area  is  attenuated  only 
about  3  db.  The  other  filters  show  comparable  results,  except  that  the 
trade-off  can  appear  as  signal  suppression  rather  than  loss  of  noise 

*Texas  Instruments,  1965,  Large  Aperture  Seismic  Array  Final 
Specifications  Rpt:  AFTAC,  Aug.,  p.  IV-42  and  IV-72. 
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rejection,  as  in  LASA  MCF-4  (Figure  VII-8).  In  this  case,  China  is  bei*'g 
attenuated  6  db  at  3.0  cps  while  the  defined  noise  areas  are  being 
attenuated  6  to  18  db.  This  attenuation  of  signal  was  observed  also 
in  MCF-1,  as  discussed  in  paragraph  C.  2.  Apparently,  the  trade-off 
under  discussion  appears  in  the  form  of  signal  attenuation  rather  than 
loss  of  noise  rejection  capability  whenever  signal  and  noise  areas  overlap 
in  £ -space,  which  is  the  case  for  both  MCF-1  and  MCF*4.  This  problem 
probably  can  be  solved  so  that  signal  is  attenuated  only  0  to  3  db  by 
specifying  a  signal -to -noise  ratio  larger  than  8.0  during  filter  synthesis, 
or  by  redesigning  the  noise  model  to  avoid  overlapping  with  the  signal  area. 

Probably  the  most  important  consideration  in  the  evaluation 
of  the  LASA  filters,  aside  from  the  signal  and  noise  k-space  response, 
are  the  filter  random  noise  responses  which  were  presented  in  Figure  VII-9. 
These  responses  indicate  that  the  output  random  noise  power  is  never 
greater  than  the  input  noise  power  between  0.5  and  4.0  cps  and,  for  all 
except  MCF-3,  input  random  noise  is  always  attenuated  greater  than 
6  db.  The  relatively  poor  random  noise  response  of  MCF-3  is  reasonable, 
as  explained  above,  on  the  basis  of  the  wide  k-space  area  coverage  of  the 
signal  model.  Since  random  noise  is  attenuated  approximately  13.98  db 
at  the  subarray  level  (assuming  all  25  channels  were  used  in  straight- 
summation  or  beam-steer  processing),  even  MCF-3  would  probably  be 
acceptable  for  processing  actual  LASA  data.  If  large  areas  such  as 
Russia  were  monitored  using  two  or  more  multichannel  filters,  each 
designed  to  cover  only  a  portion  of  the  overall  area,  a  random  noise 
response  similar  to  that  obtained  for  MCF-4  could  be  expected,  plus 
Improvement  in  organized  signal  and  noise  response  properties  would  be 
obtained. 
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SECTION  '(/III 


CPO  NOISE  AND  SIGNAL  ENSEMBLE 

A.  INTRODUCTION  ** 

This  section  discusses  the  formulation  and  processing  of  a 
noise  and  signal  ensemble  for  CPO.  The  data  used  in  this  library  was 
collected  on-site*  in  1963  in  digital  format,  using  a  DFS  (Digital  P'ield 
System)*.  This  library  of  data  is  to  be  used  in  the  following  tasks  to  be 
performed  under  this  contract: 

1)  Local  MCF  Evaluation.  Multichannel  filters  will  be  synthesized 
using  measured  noise  and  signal  for  the  respective  noise  and  signal  matrices 
This  technique  cf  filter  synthesis  will  be  compared  with  conventional 
methods  r  _'{~h  use  Lneoretical  signals  for  P  -  wave  enhancement. 

2)  Evaluation  of  Phase  Extraction.  Multichannel  filters  will  be 
synthesized  using  theoretical  coherent  noise  and  signal  models  for 
extraction  of  P,  S  and  LR  phases  of  measured  signals.  The  effectiveness 
of  phase  extraction  will  be  evaluated  'y  applying  the  MCFs  to  quarry  blasts. 

The  events  indue-  d  in  this  library  are  3  teleseisms,  12 
quarry  blasts,  1  nuclear  blast,  and  ,  noise  samples  (Table  VTII-1). 

B.  PREPARATION  OF  ENSEMBLE 

The  data  v/as  prepared  for  use  in  the  tasks  previously 
mentioned  first  by  decimation  and  then  by  whitening.  Decimation  is 
accomplished  to  yield  the  minimum  number  of  points  representative 
of  the  desiredTime -length  data  in  order  to  reduce  computer  time. 

Whitening  of  the  data  is  advantageous  for  two  reasons.  First, 
whitening  permits  maximum  use  of  computer  accuracy  in  the  wideband 
frequency  sense  when  computing  correlations  by  reducing  the  long- 
period  microseismic  energy  contribution.  Second,  whitening  is 
necessary  when  the  data  is  tc  be  used  in  time-domain  multi  .hannel 
filter  development,  since  the  ’’goodness"  of  the  filter  solution  as  a 
function  of  frequency  is  partially  dependent  upon  the  relative  spectral 
weighting.  The  steps  in  processing  this  library  are  outlined  in  the 
following  paragraphs. 

The  original  19 -channel  data,  which  is  4  min  in  length  and. 
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has  a  , 024-sec  sample  interval,  has  been  separated  into  two  groups.  The 
first  group,  decimated  by  6  to  form  144-msec  sample  interval  data, 
includes  3  teleseisms,  1  nuclear  blast  and  the  2  noise  samples  which 
preceed  and  overlap  two  of  the  teleseisms.  The  second  group,  decimated 
by  3  to  form  72 -msec  data,  includes  the  12  quarry  blasts  and  the  same  2 
noise  samples. 

Before  decimation,  each  group  of  data  was  filtered  using  a 
minimum -phase  101 -point  antialiasing  filter  appropriate  to  the  sample 
interval.  The  response  of  the  144-msec  antialiasing  filter  is  given  in 
Figure  VTII-1  and  that  of  the  72-msec  antialiasing  filter  is  in  Figure  VIII-2. 

After  decimation,  the  overlapping  noise-teleseism  pairs 
(noise  sample  AA -teleseism  BB  and  noise  sample  BB-teleseism  CC)  were 
appended  at  the  appropriate  point,  forming  two  records  of  length  470yf&~3ec . 

The  data  was  next  whitened  by  applying  an  ensemble  whitening 
(deconvolution)  filter  to  each  data  set. 

These  filters  were  designed  on  an  ensemble  basis  in  order 
that  the  relative  signal -to -noise  ratios  would  be  preserved  between  noise 
and  signal  samples,  thus  insuring  maximum  noise  and  signal  information 
which  is  necessary  in  synthesis  of  local  MCFs.  Had  each  record  been 
whitened  (deconvolved)  using  filters  designed  on  single  records,  the 
resulting  signal -to -noise  ratios  between  signal  and  noise  samples  would 
not  have  been  representative  of  the  original  data  condition. 

The  ensemble  whitening  filters  were  developed  by  computing 
the  autocorrelation  for  the  center -seismometer  output  (Z  10)  of  each  record 
in  both  data  groups  to  include  signal  and  noise.  The  autocorrelations  were 
normalized  and  stacked  (i.  e. ,  added  point  by  point)  to  form  two  average 
autocorrelations  representative  of  each  ensemble  (i.  e. ,  the  144-msec  ensemble 
and  the  72-msec  ensemble}.  Two  whitening  filters  then  were  designed 
cn  these  autocorrelations.  Each  filter  was  applied  to  its  respective  group 
of  records.  The  response  of  these  filters  is  given  in  Figures  VIII-3  and 
VIII-4.  Figure  VIII-5  shows  noise  sample  AA  before  whitening  and 
Figure  VaII-6  shows  it  after  whitening.  Figures  YIII-7  and  YTH-8  are 
before  and  after  whitening  samples,  respectively,  of  quarry  blast  CC. 

Figures  VXII-?  and  VIII-1Q  show  absolute  power-density  spectra  of  events 
AA  and  CC  before  and  after  whitening,  respectively. 


FREQUENCY  (CPS) 


Figure  VIII-6.  Noise  Sample  AA  Antialias  Filtered  and  Decimated  After  Whitening 
(Deconvolution) 
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Figure  VIH-7.  Quarry  Blast  CC  Antialias  Filtered  and  Decimated  Before  Whitening 
(Deconvolution) 


Figure  VIII-8.  Quarry  Blast  CC  Antialias  Filtered  and  Decimated  After  Whitening 
(Deconvolution) 
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Figure  VIII-9.  Power  Density  Spectra  of  Noise  Sample  AA  Before  and 
After  Whitening  (Deconvolution) 


QUARRY  BLAST  CC 
BEFORE  DECONVOLUTION 
AFTER  DECONVOLUTION 


'  A 

*  '  \  ~ 
■xV  '  V  \ 


\ 

w  -1  ' 


FREQUENCY  (CPS) 


Figure  VIII- 1G,  Power  Density  Spectra  of  Quarry  Blast  CC  Before  and 
After  Whitening  (Deconvolution) 
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SECTION  IX 


ARRAY  RESEARCH  DATA  COLLECTION 

In  support  of  TasKs  l.b.  l.b.  {2),  l.b.  (3),  and  l.d.  *  data 

have  been  acquired  from  the  Tonto  Forest  Seismological  GL-se-  /atory 
(TFO),  the  Wichita  Mountains  Seismological  Observatory  (WMO)  and 
the  Uinta  Basin  Seismological  Observatory  (UBO),  respecti  vely.  Recordings 
of  events  and  noise  samples  from  TFO  were  provided  by  AF  .'AC  with 
Texas  Instruments  assistance.  Data  from  WMO  and  UBO  were  recorded 
digitally  on-site  by  TI. 

Detai?  i  of  the  data  acquisition  will  be  provided  in  Special 

Report  No.  12. 

A.  DATA  COL]  ECTION  AT  TFO 

Data  from  TFO  for  Tasks  l.b.  (1)  and  l.bu  (2)  were 
provided  by  AFTAC  from  recordings  made  on  the  Astrodata  Digital  Data 
Acquisition  System.  These  recordings  were  made  r.i  1/2-in.  tape  hr.  an 
IBM  format  and  presently  are  being  reformatted  and  gain -corrected  at  the 
Seismic  Data  Laboratory  (SDL)  for  input  to  the  IBM  system  associated  with 
the  TIAC**  (Texas  Instruments  Automatic  Computer).  An  IBM-to-TIAC 
conversion  will  be  performed  on  this  data  to  provide  TlAC-format  data  for 
subsequent  processing. 

Responsibility  for  data  collection  at  TFO  was  given  to 
Geotech,  operator  of  the  station.  Texas  Instruments  provided  two  people 
on-site  to  coordinate  the  TFO  operation  and  to  edit  the  raw-recorded 
data  for  submission  to  f  PL.  Coordination  and  assistance  also  were 
provided  to  the  SDL  operation. 

1.  Results 

Library  data  edited  at  TFO  from  Astrodata  recordings  made 
between  16  August  1965  and  1  October  1965  comprise  the  following: 

•  Short-Period  Low  Gain  (16-19,  23-26  and  28-30  August)  ^ 
--35  moderate-to-strong  teleseisms  and  21  regional 

event  3 


*  Tasks  (modified)  as  numbered  in  contract  AF  33(657) -12747.. 
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t  Short-Period  High  Gain  (27,  JO-31  August  and 
1  September)  --24  weak  teleseisms,  35  noise 
samples,  1  regional  event  and  6  full  reels  of 
special  noise  records  (2  night-  1  early  morning, 

1  morning,  1  evening,  and  1  windy) 

e  Long-Period  Lo-’  Cain  (6,  10,  19,  20,  22,  26, 
and  29  September  and  1  October)--  8  teleseisms, 

12  noise  samples,  1  near -regional  event,  1  special 
event  and  1  weather  frontal  passage. 

®  Long-Period  High  Gain  (2,  10,  11,  17,  21,  28, 
and  29  October) — 4  weak  teleseisms,  1  regional 
event,  31  r  ,-ise  samples,  2  windy-rainy  samples, 
and  23  samples  associated  with  the  wind  and 
pressure  fluctuation  of  a  frontal  pas  save.  (Tables 
I  and  II  summarize  these  data  and  related  FDE 
information. 

It  is  reasonable  to  suspect  that  not  all  of  this  library  data 
will  be  usable,  due  to  tape  errors,  dropouts,  etc.  At  this  time,  only 
a  small  port' on  of  the  library  data  have  been  processed  through  the  Seismic 
Data  Laboratory  (SDL),  and  that  has  not  been  evaluated  yet. 

A  standard  8 -min  edit  length  was  used  f^r  all  short -period 
data.  Edits  from  full  reels  of  special  noise  tapes  were  made  to  include 
a  30 -sec  overlap,  should  reconstruction  of  an  entire  reel-length  sample  be 
desired  later.  Long -period  edits  were  40  min.  Since  field  reels  could 
run  from  65  to  7°  min,  it  was  possible  to  obtain  two  40-min  samples  from 
one  reel  with  a  10-  to  15 -min  overlap. 

Pertinent  daily  calibrations  were  edited  for  submission  to 
SDL  for  seismometer  channel -gain  correction.  These  calibrations  as  well 
the  short-  and  long-period  irequency  responses  are  to  be  analysed  and 
the  results  are  to  be  made  available  t  TI. 

Included  in  the  librory  for  the  purpose  of  quality  control  are 
records  of  d-c  pulses  for  long-  and  short-period  channel  verification,  and 
records  of  system -noise  tests  made  with  dummy-loade  i  PTA*  and  shorted- 
out  Voltage -Controlled  Oscillator  (VCO)  inputs**  at  the  remote  vans. 


*  Dummy-loaded  PTA  testing  involved  replacing  the  seismometer  by  a 
resistor. 

**  Shorted -out  VCO  testing  involved  removing  PTA  fvorn  telemetering 
circuit  and  shorting  across  the  VCO  input. 


SUMMARY  OF  FDE  INFORMATION 


cO  0s 

r- 

m  —i 

m 

o 

(M 

(M 

fn 

o 

00 

CO 

pH 

CO 

CO 

O'*  ~h  rp 

oj 

o  r- 

o 

cc 

00 

rc 

cc 

O' 

s£> 

(M 

O 

00 

O'  ^ 

O' 

o- 

"1  rl  S'  O' 


NhOo'ODt-rtfnO'O'-i  i  O  ifl 


m  to  tc  rt*  m  i  m  m 


m  (O  M  N 
Cl  «  rt 


PO 

O'  o  ^ 

H  H  III 


w  ps  a?  # 


OOO  fOOMHOOIJ'ffllfl 

—I 


50  05 

U  m  +j 

*2  ~  3 

£  O  0 

f,  fll  0  I 

u  ^  o 
ffi  tj  «  U 

*  C*  -fi 

.«  g  e  « 


w  w  f> 

00  N  N 


o  a)  o  £! 
fQ  2  &  0 


W  £  £ 

vC  *-H  lO 


?  -s 

«  3  M 

m  oi  rQ 

o  a  o> 

fi  ¥  IB 

43  <u  *  £ 
**+>►. 

M  cd  S  ;£. 
O  S  O  .? 

Z  OZh 


W  05  0J 

05  05  05 

-a  T)  T3 


+1  35  0) 

£>  Z  Z 


tO^OnfOfOOvOsOf-NvONNd1 


vO  in  (\j 
o  O'  r-  O' 


CO  £  £  2 
m  n  c  n 

•  ti» 

r-  in  rj  in 

r¥  if)  — I 


O'*  -*  t- 

*  •  i  -fr 

N  ^  N  O' 
05  O  —<  ' — i 
C"  in  -^3  (M 
lO  (O  H  o 
NOW'? 
nj  i-¥  — i  f- i 


voa-'c-  co(*iinininMOOt"OMOi»N 

NO  ^  p~  vfi  t»  ?•  N  f-  vO  ?•  >0  t>  ^  ¥  *o  >o 


in  os  Z  tn  tn  tn  in  ui  w  Z  Z  in  m  tn  Z  in  <’i 

eg  — i  n  '?o'ptrot>owo'Nin?-'!,(j'f' 

in— iuri  ot^mmrnvo— 'pnoo'inr^vOvo 

¥  N  NWNNN-HOOHrtrfrftOHH 


in  h-o  in  ¥  so  n  in  m  t*-  on  rj 

•  ••  ******^Sr-K'  »  t  , 

h  n  ,n  Nffl®oO'0'CPd|¥ininwmN 

ini?NNNNNinoNN¥oo 
r-  ifi  in  oo¥i?ini3j¥inNinin¥ro>oN 

.— i-^OM  — (  ©  — iNWinNWWNWiCUIW 

■O  O  H  NO^'?i?'?inNNr-(OOOrO 

— *  no  nj  nj  — *  — *  ■ — ■  — *  - — ■  > — t  ©  o  — *  — *  o  ©  o 


WibOtfltC  bfl  Wl  <50  bCbObO&ObObObOtiOSiebOhOlaOSflbfi 


<*,<<<  <<<  <<<<<<<<<<<;<<:<: 

'Or-r^n-  s-  c-  r-  Ntooseocofflcoo'o'o'o'ooo 

*“t  H  »H  »H  »H  »— «t  M  (H  »“<  *H  »H  M  »*—f  (VJ  f\J  fVj  N  CO  CO  CO 


Short-Period  High  Gaia 


cO 

m 

fH 

»-H 

O  n£> 

p 

CM 

fM 

CM 

CO 

»—<  Is- 

O' 

co 

vO 

CO 

vO 

fH 

• 

• 

• 

• 

* 

in 

m 

P 

•4* 

in 

Ph 

£ 

co 

CM 

CO 

O 

04  cn 

00 

O' 

O' 

CO 

co 

r- 

cO 

cO 

CM  CO 

cO 

CO 

cO 

CO 

^H 

P 

i4  c|  “  Q.-P  u 

bo  .-i  nj  sn  nl 

C  C  «  .2  nj  X 

O  O  O  &  o  nJ 


03 

1-1 

p 

o 

03 

IH 

at 

at 

4-> 

at 

at 

<y 

«h 

CO 

r-H 

0) 

CO 

03 

VH 

0 

0, 

<\i 

<4-t 

O 

G 

tH 

Pi 

0 

03 

CO 

►-H 

hH 

’(-< 

•M 

•  J 

<1) 

i 

o 

u 

M 

43 

IH 

o 

> 

3 

O 

a 

M 

(D 

r* 

S* 

03 

O 

TJ 

•H 

c 

P"t 

V 

0) 

#H 

ni 

0t 

I-H 

at 

X* 

03 

n: 

cn 

at 

0 

ft 

pH 

r-* 

43 

03 

l—l 

nJ 

0} 

at 

0 

g 

< 

03 

1— 1 

n) 

•rt 

T) 

-M 

O 

+J 

X 

p* 

at 

Jh 

•H 

■  r* 

O 

0) 

be 

f  ^ 

4-> 

bo 

o 

Vi 

O 

9} 

<D 

3 

c$ 

CO 

at 

cv 

V 

3 

pM 

o 

0< 

at 

0 

o 

Pi 

at 

N 

d 

£ 

03 

a 

3 

O 

3 

o 

<+H 

IM 

<3 

a) 

3 

O 

2 

n 

£ 

> 

w 

Vh 

ffi 

►o 

tn 

S 

z 

o 

£ 

h 

0 

O 

H 

£  U  £ 

£ 

£ 

£ 

U 

w  w 

W 

w  w  w 

£  £ 

£ 

N  vO  N 

O 

00 

in 

m 

O'  rM 

CO 

CM  P  O' 

O'  P 

On 

•  •  • 

• 

• 

• 

♦ 

•  • 

• 

•  •  • 

•  * 

« 

\1  N  M 

vO 

CM 

On 

00 

O  P 

m  t--  oo 

C4  P 

M’ 

N  't  O' 

O' 

O' 

CM 

cO 

p  p 

P  'O  m 

t—  CO 

r- 

£ 

Z  Z 

z 

Z 

w  £ 

co  £ 

CO 

00 

co  in 

o 

m  <m 

in  vO 

fH 

• 

On 

CO 

•  • 

O'  co 
co  p 

<M 

iH 

• 

p- 

*— H 

«  • 

f-3  CO 

•  • 

00  NO 

r— 1 

• 

'M 

CM 

p 

cO 

in 

CM 

O' 

*55- 

in 

m 

v 

Jt 

# 

*• 

w 

• 

X 

• 

• 

• 

• 

• 

t— 

CO 

m 

P 

H 

00 

b- 

n- 

ah 

p 

n£> 

o 

o 

CM 

cO 

m 

CM 

o 

P 

CM 

o 

P 

m 

CM 

fH 

cO 

CM 

•O 

cO 

m 

rO 

o 

fH 

m 

CM 

O' 

1— 1 

O'- 

00 

*-4 

cO 

00 

CM 

cO 

cO 

'O 

cO 

p^ 

in 

cO 

P 

o 

CO 

CM 

m 

p 

cO 

O 

cO 

fH 

CM 

NO 

t- 

fN- 

00 

00 

fH 

fM 

b- 

r- 

On 

O' 

n£> 

b- 

r- 

^H 

i—i 

o 

o 

o 

o 

o 

o 

O 

o 

o 

O 

•■“I 

O 

o 

o 

CM 

O 

so  bo  bo 
3  3  0 
<  <  < 

t*-  f'*  ts- 

rM  fM  fM 


*  -*->  4-> 

4->  -M 

5®  a  cv 

cu  a 

a 

P  03  iy 

V  <D 

<u 

<J  cn  co 

in  w 

CO 

f-H  p-H  p-H 

f— »  sD 

O' 

CO 

H 

IX -5 


ng -Period  Low  Gain 


o 

u 

m  I  ° 

T?  4)  {I 

°  2  m 

4->  ^ 

CO  «“* 

*  id  £ 

O  y  ^ 


c 

O  g. 

+->  nJ 

os  1-71  rt 

rt  c  ‘3 
O  Q  3  j!  £ 

u  u  x  J!  c 


O  (tf  Jn-ri+j  m  ij  jjj 
u.  ^  ®  fi  P  H 

wn)i)®nJOn(j 

oo^Swmwu 


oo  o  n 

H  Tf  00 


^  T}<  ^ 


The  desired  quantities  of  various  types  of  data  were  not 
obtained  always.  The  effective  recording  period  was  shorter  than  anticipated. 
Few  long  “period  events  that  had  a  P,S  and  L  signal  -to-noise  ratio  greater 
than  1  which  did  not  over -modulate  the  recording  system,  vere  recorded. 

No  regional  events  with  known  epicenter  were  observed,  since  station 
personnel  had  not  determined  areas  of  mining  activity  greater  than  aear- 
regional  distances.  There  was  no  frontal  passage  during  short -period 
recording,  nor  was  there  any  rain  not  accompanied  by  lightning  and  spiking. 

B.  DATA  COLLECTION  AT  WMO 

Data  collection  in  support  of  Task  l.b.  (3)  was  begun  at  the 
Wichita  Mountains  Seismological  Observatory  (WMO)  on  16  November  1965. 
Data  were  recorded  from  3 -component  strain  instrumentation,  3-component 
inertial  instrumentation  and  the  13 -channel  surface  areal  array. 

It  should  be  noted  that  the  cutoff  date  of  this  report  precluded 
inclusion  of  the  details  of  WMO  data  collection,  however,  the  task  was 
completed  at  the  time  this  report  was  published. 

A  later  report  will  detail  the  results  of  the  data  recording 
and  wil)  present  a  convent' onal  seismological  analysis  of  recorded  events. 
Noise  samples,  events  and  special  tests  selected  by  the  analyst  for  inclusion 
in  the  library  will  be  transcribed  on  IBM -format,  tape  and  will  be  supplied  to 
the  Vela  Seismic  Center, 

It  was  planned  to  record  for  a  total  of  nine  days  during  the 
period  ending  30  November.  In  general,  normal  station-working  hours 
were  to  be  followed  especially  in  view  of  the  remote  location  of  the  station 
and  the  accompanying  military  activities. 

It  was  planned  to  record  the  following  data: 

•  Twenty  30 -min  ambient  noise  samples 

•  10  strong  teleseisms 

«  20-30  quarry  blasts  and/or-  near-regionals  (with  well- 
developed  surface  wavetrains) 

e  Calibrations  (including  all  daily  and  special  calibrations  ) 

•  System-noise  tests  (at  start  and  finish  of  data  collection) 

All  data  were  to  be  field-edited.  The  data  on  each  reel  were 
to  be  preceded  by  a  brief  DFS  calibration  and  5-min  recording  (truck- 
noise  test  made  with  the  operational -amplifier  inputs  grounded).  Skew 
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tests  were  to  be  performed  at  least  once  per  day. 

The  following  DFS  channelization  was  to  be  used; 
Channel  Contents 


2 

3 

4 

5 

6 

7 

8 
o 

10 
11 
12 
13 -2  b 
26 

27 

28 

29 

30 

31 


Vertical  Strain--SZS--High  Gain 
Horizontal  Strain --SNS- -High  Gain 
Horizontal  Strain --SES- -High  Gain 
Vertical  Inertial --SPZ --High  Gain 
Horizontal  Inertial--SPN--High  Gain 
Horizontal  Inertial- -SPE --High  Gain 
Vertical  Strain--SZS  -  -Low  Gain 
Horizontal  Strain --SNS- -Low  Gain 
Horizontal  Strain--SES--Low  Gain 
Vertical  Inertial --SPZ --Low  Gain 
Horizontal  Inertial--SPN--Low  Gain 
Horizontal  Inertial- -SPE --Low  Gain 

13 -Element  Surface  Array - Low  Gain 

SPN  at  Z6  Location 

SPE  at  Z6  Location 

Summation,  13  Surface  Verticals 

Spare 

WWV 

Station  Time 


Channels  1-6  were  to  be  recorded  at  such  gain  settings  that 
would  cause  typical  noise  peaks  to  come  within  10  db  of  DFS  clipping.  Each 
channel  was  to  be  adjusted  individually  in  order  to  make  maximum  use 
of  dynamic  range.  Channels  7-12  were  to  have  gains  in  the  same  proportion 
as  channels  1-6  ,  but  approximately  20  db  lower.  The  amount  of  decreased 
gain  was  to  be  determined  by  whatever  was  required  to  record  reasonably 
large  teleseisms  without  digital  clipping.  Channels  13-25  were  to  be 
operated  with  a  gain  sufficiently  low  that  average  teleseisms  could  be 
recorded  without  clipping. 
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Station  develocorder  film  and  all  logs  were  to  be  obtained 
from  the  station  manager  and  sent  to  Dallas  for  copying.  Film  and  logs 
from  the  strain-instrumentation  group  were  to  be  obtained  from  Garland. 

The  station  manager  and  the  strain-instrumentation  group 
leader  were  to  be  kept  informed  of  TI  on-site  activities.  Efforts  were  to 
be  made  to  minimize  departures  from  station  routine, 

C.  DATA  COLLECTION  AT  UBO 

1.  Program 

Task  l.d.  of  the  program  requires  that  a  Seismic  Data 
Library  be  created  from  on-site  digitally  recorded  data  collected  at  the 
Uinta  Basin  Seismological  Observatory  (UBO).  These  field  data  are  to 
be  edited,  gain-corrected  and  made  available  to  AFTAC  on  IBM  tapes, 
together  with  a  report  of  a  conventional  seismological  analysis. 

2,  History 

The  AFTAC  Digital  Recording  System  was  dispatched  from 
Dallas  to  UBO  on  13  September  1965,  but  accidentally  was  overturned  at 
Moab,  Utah,  on  16  September,  This  truck  eventually  was  repaired  enough 
to  be  driven  back  to  Dallas.  Because  the  program  at  UBO  involved  collecting 
data  which  would  not  be  available  at  a  lat»r  date,  Texas  Instruments 
made  ready  and  dispatched  its  own  Digital  Recording  System  to  UBO  to 
accomplish  the  date,  collection.  This  second  truck  departed  Dallas  on 
20  September  and  arrived  at  UBO  on  22  September.  All  data  sources  being 
recorded  w’re  cabled  into  the  Digital  Recording  System,  and  internal 
tests  and  calibrations  were  performed.  The  first  recordings  with  all  input 
systems  properly  calibrated  and  operating  were  made  on  3  October  1965. 

Several  data  sources  had  to  be  operational  simultaneously 
for  the  ultimate  schedule  of  instrument  evaluation,  array  evaluation  and 
eventual  library  formation  to  be  successful.  Sources  required  were  the 
two  short -period  areal  arrays,  both  multichannel  filters,  and  the  6 -element 
vertical  array.  The  latter  array  was  a  temporary  installation  whose 
operation  required  that  all  other  schedules  be  keyed  to  its  availability. 
Consequently,  all  time  and  equipment  schedules  were  carefully  coordinated 
among  Texas  Instruments,  AFTAC  and  Geotech. 

Recording  continued  from  3  October  on  a  routine  schedule 
until  the  6 -element  vertical  array  malfunctioned  on  S  October.  Investigation 
showed  extensive  and  time-consuming  repairs  would  be  necessary  before 
usable  data  would  be  available  again.  This  time  lag,  other  commitments  of 
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the  Texas  Instruments  Digital  Recording  System,  and  the  accumulation 
of  a  satisfactory  amount  of  ..sable  recorded  data  from  UBO  prompted 
the  decision  to  return  che  system  to  Dallas  on  14  October, 

3.  Data  Sources 

Data  sources  and  their  channelization  in  the  Digital 
Recording  System  are  given  in  Table  I.  Outputs  from  the  multichannel 
filters  were  taken  parallel  to  the  Develocorder  input  terminals.  Outputs 
from  all  UBO  instruments  were  obtained  at  the  PTA  outputs;  outputs  from 
the  6-element  vertical  array  were  taken  from  the  Signal  Conditioner  and 
Filter  package  located  in  the  auxilliary  instrument  van.  During  all  sine- 
wave  calibrations,  a  signal  was  obtained  directly  from  the  driving  function 
generator  and  recorded  in  channel  30. 

4.  Recorded  Data 

The  recordings  at  UBO  wer;  made  on  17  field  tapes  dis¬ 
tributed  as  shown  in  the  following  list: 

©  High-gain  noise  6  reels 

«  Normal-gain  events  8  reels 

•  Calibrations  2  reels 

©  Tests  1  reel 

The  six  high -gain  tapes  contain  a  total  of  11  hr  40  min  of 
noise  data  recorded  in  12  separate  sessions.  The  eight  event  tapes  contain, 
on  the  bas5s  of  a  preliminary  analysis,  15  lecords  of  either  quarry  blasts 
or  near -regional  earthquakes,  and  25  teleseisms. 

Calibration  and  test  reels  contain  daily  calibrations  and  special 
frequency -response  calibrations  and  system-noise  tests. 

5.  Utilization  of  Data 

Approximately  half  of  the  data  recorded  at  UBO  have  been 
edited  and  transcribed  into  TIAC  records.  This  edit  was  necessary  to 
provide  data  for  an  early  evaluation  of  the  vertical  array.  Details  of  this 
program  are  given  elsewhere  in  this  report. 

The  remaining  data  are  being  edit  ,-d  currently  and  will  be 
transcribed  into  TIAC  records.  An  analysis  of  array  calibration  data  will 
/ield  the  proper  multiplication  factors  to  normalize  each  data  channel  to 
some  uniform  equivalent  transducer  gai-\  This  gain-corrected  TIAC 


ensemble  will  be  transcribed  next  into  an  IBM -tape  library  which  will 
be  made  available  to  AFTAC.  A  listing  of  library  contents,  along  with 
a  conventional  seismological  analysis,  will  be  presented  in  a  special 
eport. 
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SECTION  X 


DISCUSSION  OF  LASA  PROCESSING  REQUIREMENTS 

A.  INTRODUCTION 

At  the  invitation  of  the  Advanced  Research  Projects  Agency 
in  September,  Texas  Instruments  personnel  prepared  and  delivered  several 
presentations  at  the  LASA  Coi  ference  in  Washington,  D.C.  One  of  these 
presentations  dealt  with  the  problems  associated  with  processing  LASA 
data.  A  brief  summary  is  given  here. 

The  full  text  of  this  paper  will  be  published  as  Special  Report 

No.  13. 

B.  GENERAL  CONCEPT  AND  PROCESSING  FUNCTIONS 

1,  On-Line  Detection  and  Approximate  Location 

The  total  processing  of  data  from  a  LASA  and  the  nature 
and  quantity  of  data  permanently  stored  on  magnetic  tape  are  determined 
by  the  status  of  the  seismic -wave  field  at  the  time.  For  example,  when 
ambient  noise  is  present  with  no  teleseismic  events  above  a  given  magni¬ 
tude  threshold  and  the  absence  of  such  events  is  known,  no  further 
processing  is  required  and  only  a  limited  amount  of  data  need  be  stored. 

In  the  presence  of  a  teleseismic  event  above  the  threshold,  additional  on-line 
processing  dependent  on  the  approximate  location  of  the  event  would  be 
carried  out,  A  greater  amount  of  data  would  be  stored  in  the  presence 
of  a  teleseismic  signal  with  certain  characteristics.  In  the  presence  of  an 
abnormally  high  noise  field,  a  large  amount  of  the  data  (possibly  up  to  the 
525  channels)  would  be  stored  so  that  maximum  processing  capability 
might  be  exploited.  Therefore,  the  major  on-line  function  ol  the  processor 
must  be  determination  of  the  status  of  the  seismic -wave  field.  The 
processor  must  detect  teleseisms  above  a  specified  magnitude  threshold, 
but  need  determine  their  location  only  to  sufficient  precision  to  permit 
additional  on-line  processing  for  extracting  more  detailed  information  about 
the  event. 

2,  Data  Storage 

An  on-line  temporary  buffer  storage  for  the  full  525  channels 
is  contemplated.  The  buffer  storage  must  save  the  raw  data  from  the  525 


channels  long  enough  for  appropriate  action  to  be  taken  fo^*  detected  events. 
This  period  is  probably  on  the  order  of  1-1/2  min. 


Magnetic  tape  units  are  provided  for  permanent  stora^h. 
Worldwide  P-wave  outputs  from  the  21  subarray*  would  be  stored  for  all 
data.  On  the  basis  of  the  detection  logic,  a  number  of  other  storage 
operations  up  to  and  including  storage  of  the  full  525  channels  would  be 
provided. 

3.  Postdetection  Processing 

Postdetection  Processing  would  include  immediate  on-line 
processing  of  the  buffered  data  to  develop  the  best  steered  beam  for  each  of 
the  21  subarrays  based  on  the  approximate  location  of  the  event.  Additional 
processing  of  this  data  would  include  improved  location  and  measuiement 
of  arrival  time  and  would  accomplish  signal  extraction  for  subsequent, 
measurement  of  classification  parameters.  In  the  case  of  larger  teleseisms, 
special  redetection  processing  would  be  provided  for  locating  possible 
smaller  events  which  might  be  buried  in  a  larger  teleseismic  arrival. 

4.  Built-in  Learning  Function 

It  is  anticipated  that  determining  the  necessary  corrections 
(particularly  between  subarrays)  to  provide  maximu..:  enhancement  in  the 
beam -f c  rming  function  will  require  a  long  period  of  time  and  the  use  of 
data  from  a  large  number  of  teleseisms.  Large -detected  events,  therefore, 
would  undergo  processing  for  upgrading  of  the  beam-  former  corrections. 

5.  Off-Line  Network  Processing 

Finally,  the  function  of  combining  information  from  a.  number 
cf  LASA  stations  must  be  performed  tc  accomplish  source  classification. 

6.  General  Implementation 

A  generalize  i  picture  of  this  processing  concept  is  shown  in 
Figure  X-l.  The  325  channels  from  the  subarrays  are  delivered  to  a 
special-purpose  processor  which  provides  signal  detection  and  preliminary 
location.  In  addition,  the  special-purpose  processor  provides  21  processed 
broadband  subari*ay  outputs  for  permanent  magnetic -tape  storage.  The 
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SECTION  XI 


AUTOMATED  MAPPING  SYSTEM 

A.  INTRODUCTION 

A  series  of  five  programs  has  been  written  for  the  IBM- 
7044  Computer  and  the  Cslcomp  Plotter  to  aid  in  the  evaluation  of  potential 
array  sites  and  in  the  interpretation  of  data  recorded  at  an  array  site. 

These  programs  utilize  sets  of  worldwide  geographic  coordinates  and  array - 
site  coordinates,  and  provide  maps  of  various  types  including  stereographic 
projections  centered  on  the  array  site  and  a  map  of  the  world  in  k-space 
centered  on  the  site. 

Descriptions  of  these  five  programs  and  information  required 
to  run  each  program  will  be  published,  together  with  a  listing  of  Fortran 
statements,  as  Special  Report  No.  11.  Source  decks  for  these  programs 
and  for  the  Calcomp  Plotter  subroutines  have  been  supplied  to  the  Project 
McniL  at  the  Vela  Seismic  Center. 

To  obtain  coordinate  seta  of  the  major  land  masses,  seismic 
areas  and  other  points  of  interest,  a  Polar  Stereographic  Projection  map, 
published  by  the  Aeronautical  Chart  and  Information  Service,  U.  S.  Air 
Force,  having  a  scale  of  a  3roxirr.ately  1:25X10^  was  used. 

Approx:mately  2900  points  were  used  to  outline  continental 
and  political  bo”ndaries  and  major  islands.  An  additional  460  points  were 
required  to  show  major  seismic  areas;  266  points  were  used  to  show  .5° 
increments  of  latitude  and  longitude..  About  3600  points  are  contained  in  the 
present  library. 

Cards  and  listings  containing  these  p  in  have  been  supplied 
to  the  Project  Monitor. 

B.  APPLICATION  OF  MAPPING  SYSTEM 

This  mapping  system  provides  an  economical,  graphical 
solution  to  the  following  problem:  from  a  specified  site,  what  is  the 
appearance  of  the  world  in  k-space  (for  a  specified  frequency  and  propagation 
phase)  nd  in  geographic  space?  In  general,,  the  programs  will  plot  seta 
of  coordinates  as  any  definable  function  of  distance  from  any  one  site.  By 
plotting  major  e<n Tnquake  epicenter  areas,  it  'an  be  determined  if  a  site 
is  likely  to  record  a  large  number  of  earthquakes  or  If,  perhaps,  it  li^s 
in  a  shadow  zone  for  many  seismic  area  By  jotting  possible  areas  of 
nuclear  activity,  it  can  be  determined  if  a  site  will  .ikely  record  events 
from  a  majority  of  them.  By  plotting  coordinates  of  the  reported  eaith- 
quakes  of  tne  last  few  years,  one  can  determine  if  there  is  an  "earthquake 


Xl-i 


screen"  between  the  site  and  a  possibly  area  of  interest. 

For  use  in  evaluation  or  interpretation  of  several  sites,  a 
graphical  solution  may  be  obtained  by  plotting  the  same  data  centered 
on  each  site  and  comparing  individual  results  simultanously.  Plans  are 
being  formulated  now  for  representing  data  from  several  sites  on  the 
same  chart.  This  capability  will  be  valuable  for  multiple -site  evaluation 
and  interpretation. 

1.  Examples 

Figures  XI -1  through  XI -3  show  the  finished  output  from  the 
series  of  programs.  Data  were  taken  from  the  Calco.np  Plotter  and  the 
maps  were  drawn  by  a  draftsman  using  the  plotted  ooints. 

Figures  XI-1,  XI -2  and  XI -3  serve  a  3  fold  purpose. 

Figures  XI-1  contains  the  following  information: 

•  It  is  a  stereographic  projection  of  the  areas  of  the  world 
90°  or  less  from  Tonto  Forest  Seismological  Observatory 
(TFO),  near  Payson,  Arizona.  True  distance  and  azimuth 
may  be  measured  from  TFO  to  any  point  on  the  nr  up, 

®  It  presents  15°  latitude  avid  longitude  line 3  (shown  as 

screened  light  grey  lines)  as  they  appear  on  the  sterographic 
projection 

«  It  shows  the  earth's  major  seismic  areas  (black  triangles) 
as  seen  on  the  steroographic  projection  from  TFO. 

Figure  XI. -2  shows  the  following  information: 

®  It  is  a  storeographic  projection  centered  on  the  antipode 
of  TFO  or  that  point  diametrically  opposite  TFO;  this 
encompasses  all  points  greater  than  90°  from  TFO. 

True  distance  from  TFO  to  any  point  on  this  map  may  be 
computed  by  finding  the  distance  from  the  antipode  of 
TFO  to  the  print  and  subtracting  it  from  1  80°.  Azimuth 
from  TFO  to  any  point  on  this  map  can  be  figured  by 
finding  the  azimuth  from  the  antipode  and  subtracting  it 
from  360°. 

•  It  presents  15°  latitude  and  longitude  lines  (shown  as 
screened  light  grey  lines)  as  seen  on  the  stereographic 
projection  centered  on  I/O's  antipode. 
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•  It  shows  the  earth's  major  seismic  areas  (black  triangles) 
greater  than  90°  from  TFO,  as  seen  on  the  stereographic 
projection. 

Figure  Xl-3  presents  the  following  information: 

9  It  is  a  map  of  the  world  as  a  function  of  wavenumber 
seen  from  TFO  for  a  1 -cps  P-wave.  The  points  plotted 
lie  on  true  azimuths  from  TFO.  The  shadow  zone  contains 
all  points  greater  than  104°  from  the  station,  that  being 
the  approximate  distance  beyond  which  initial  P-wave? 
are  not  recorded.  Points  actually  located  near  TFO  are 
greatly  distorted  on  the  projection  because  of  their  large 
k-number;  thus,  the  map  shows  only  points  of  Ic  =  .  09  or 
less. 

•  It  presents  15°  latitude  and  longitude  lines  (shown  as 
screened  light  grey  lines)  as  seen  on  the  It -projection 
from  TFO. 

•  It  shows  the  earth's  major  seismic  areas  (black  triangles) 
in  £-  space. 
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Figure  XI-3.  World  Map  Seen  fr«.m  TFO  for  a  1.0  CPS  P-Wave  Showing 

K-Proiection  Lat.  and  Long.  Lines,  and  Earth’s  Major  Sei3rnic 
Areas  m  K-Space 


